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ABSTRACT 


A new type of pressure transducer, pressure-sensitive paint, was used to obtain 
pressure distributions associated with shock-boundary layer interaction. Based on 
the principle of photoluminescence and the process of oxygen quenching, pressure- 
sensitive paint provides a continuous mapping of a pressure field over a surface of 
interest. The data measurement and acquisition system developed for use with the 
photoluminescence sensor was evaluated first using an underexpanded jet blowing 
over a flat plate. Once satisfactory results were obtained, the system was used to 
examine shock-boundary layer interaction in a blow-down supersonic wind tunnel 
at Mach numbers of 1.4 and 1.7. Details of the measurement technique, and 


discussion of the flow fields which were examined, are reported. 
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I. INTRODUCTION 


One of the major thrusts of military aircraft engine design is towards achieving higher 
pressure ratios per compressor stage in order to reduce the number of stages and thereby 
reduce engine weight. This requires higher aerodynamic loading on the blades, and/or 
higher relative flow velocities over the blading. Where the relative Mach number is 
supersonic, shock waves will occur ın front of and within the blade passages. The 
interaction of strong shocks with blade surface boundary layers can lead to high losses and 
reduced turning. Thus shock-boundary layer interaction, and methods to alleviate its 
attendant negative effects on compressor blade performance, warrant careful investigation. 
It was the need to make surface pressure measurements in a model simulation of shock- 
boundary layer interaction in a fan blade passage, which motivated the present study. 

Currently, aerodynamic (surface) pressure measurements are made by installing a 
prearranged distribution of pressure ports in a wind-tunnel model or prototype aircraft. 
Since the installation of taps is an expensive process, the location and number of taps 1s 
normally predicated on a pre-selected area of interest. Obviously, this arrangement offers 
little flexibility if the area of aerodynamic interest is broad and fine spatial resolution is 
needed. Even with a substantially instrumented region, the information obtained is 
spatially discrete. It is often difficult to decide where to position the taps before the model 


is tested. 


Within the last six years, the Central Aero-Hydrodynamic Institute in Moscow, the 
University of Washington (with support from Boeing and NASA) and McDonnell Douglas 
Aerospace have independently developed a new technique to measure surface pressures 
which is based on photoluminescent coatings [Ref. 1 through Ref. 9]. The pressure- 
sensitive coatings contain photoluminescent probe molecules. With the proper excitation 
of light, the molecules luminesce with an intensity which is inversely proportional to the 
local surface pressure. As a result, pressure-sensitive paints offer almost limitless spatial 
resolution, in addition to being applicable to surface locations which may not be accessible 
for pressure ports. 

A. DESCRIPTION OF THE CONCEPT 

Figure 1 illustrates the general concept of the method. A luminescent molecule 
dissolved in an oxygen-permeable plastic resin is applied to a surface of interest. When 
illuminated by ultraviolet radiation the coating luminesces with an intensity that depends 
on the oxygen partial pressure seen by the active molecule. Thus the oxygen partial 
pressure can be determined from the emitted light intensity and since the mole fraction of 
oxygen in air is a known constant, the air pressure may then be readily calculated. The 
result of this oxygen dependence 1s that the airflow-induced surface-pressure field (and the 
corresponding oxygen partial-pressure field), gives rise to a luminescence intensity 
distribution. The intensity field can be imaged using a video camera, the image digitized, 


processed and stored on a computer. The end result is a map of the surface pressure field. 
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Figure 1. Schematic Representation of the Pressure-Sensitive Luminescent Coating 
Concept! . 


The processing of the surface pressure image relies on a calibration curve for the coating 
material, which 15 based on the Stern-Volmer relation.[Ref. 1: p. 34] 
B. APPROACH TO THE PROBLEM 

The present work was focused on the development of a working 
pressure-measurement system, utilizing pressure-sensitive paints, for use in the Gas 
Dynamics and Turbopropulsion Laboratories at the Naval Postgraduate School. Although 
the system's initial application was tailored for use in the study of shock-boundary layer 
interaction in a pilot transonic cascade wind tunnel, and subsequently in a larger 


supersonic wind-tunnel, its use for other compressible flow studies was anticipated. 


! See Ref. 1, page 34, Figure 1 


The initial development was carned out using an existing flat-plate model in an 
underexpanded sonic jet. This was an arrangement that would produce effects similar to 
those to be measured in the tunnel, (i.e. shock waves impinging on a surface, creating a 
near-discontinuity in a pressure distribution) but be visually unrestricted since no walls and 
windows were present. Subsequently, measurements were made of the interaction of the 
starting shock wave with the side-wall boundary layer in the test section of a blow down 
wind tunnel. Results were obtained at M=1.4 and M=1.7 using different fixed nozzle 
blocks. The blocks for M=1.7 were designed and built for the present study, and the 
design is detailed in Appendix E. 

In reporting the work, Section II first provides some preliminary background work 
done in photoluminescent barometry in addition to the theory of photoluminescent 
coatings. Section III describes the development program using the sonic free-jet. Section 
IV provides results of applying the system in the supersonic wind tunnel, and Section V 
gives conclusions and recommendations. Since the system represents a powerful new tool 
for two laboratories, the attempt has been made to include all information and details 


necessary to both use and extend the present system. 


II. THEORY OF MEASUREMENT 


A. BACKGROUND 

Photoluminescence and the process of quenching by oxygen is a photo-chemical 
interaction that is well documented and reasonably well-understood. In 1980, Peterson 
and Fitzgerald first recognized its potential as a means to make qualitative assessments in 
flow visualization studies [Ref. 3]. They reported on a technique in which fluorescent dye, 
when applied to a surface and excited by a blue light, illuminated in varying degrees of 
intensity, depending on the introduction of nitrogen or oxygen. When nitrogen was 
injected through a wall static-pressure tap, a bright streak of luminescence, moving in the 
direction of the surface flow, was observed. The nitrogen, in close proximity to the 
fluorescent dye, reduced the effect of the quenching (by displacing the oxygen) and 
allowed the intensity of the luminescence to increase. The introduction of oxygen, on the 
other hand, enhanced the quenching process and decreased the luminescence, prodieimg a 
dark streak. This demonstrated that the concept of quenching could be used to acquire 


surface pressure measurements. 


B. LUMINESCENCE 


According to the laws of thermodynamics, a system in an excited energy state due to 
the absorption of radiation, must be such that the reverse process of emission of radiation 


can occur. Ifno other process is available to allow the system to return to its initial state, 


then emission must take place. Light emission excited by light absorption is called 
photoluminescence. 

Luminescence is a broad term which encompasses both fluorescence and 
phosphorescence. The basic, underlying difference between these two processes is the 
time of emission with respect to the end of excitation; fluorescence has a relatively short 
lifetime, approximately 107-10" seconds, whereas phosphoresence's lifetime persists for 
107-107 seconds. It is the phosphorescent emission that was used in the present work to 


quantitatively measure the static pressure on a surface of interest. 


C. ABSORPTION AND EMISSION MECHANICS 

Quantum mechanics provides the necessary theoretical tools to describe the absorption 
and emission processes of light that were previously unexplainable by wave theory 
[Ref. 10: p. 457]. Instead of considering light as a wave train of infinite length, the 
radiation is instead composed of discrete photons (or quanta) whose energies are equal to 
hv, where h is the universal Planck constant, 6.63 x 10” erg sec. The concepts of light 
being consider as a wave train and that of discrete photons are linked in that the energy of 
the quantum is proportional to the light excitation frequency, v . This led to the 
hypothesis that radiation of frequency v could not be emitted or absorbed in arbitrary 
amounts, but only in quanta of energy € — hv. During an emission process following light 


absorption, less energy release by emission could occur, but would radiate at a lower 


frequency. 


An absorbing molecule can exist in discrete, sharply defined levels in which transitions 
between states occur in certain definite steps. The lowest of these states is the ground 
state. States of higher energy are called excited states. There is a large number of 
possible excited states, but usually only several ofthose of the lowest energy are produced 
with visible or ultraviolet light. A simplistic description of the nature of these excited 
states is that a single electron has been promoted from an orbital occupied in the ground 
state to an orbital that is vacant while in the ground state. 

Figure 2 illustrates several energy levels of an arbitrary molecule which are represented 
by the horizontal lines. The vertical distance between two of these lines is proportional to 
the corresponding difference in energy, hv; the level N represents the ground state. By 
absorption of light frequency Vnr the molecule is raised to the state F and if no other 
energy levels exist between N and F, the molecule can return to N only by re-emission of a 
photon of light ofthe same frequency. Theoretically this is the simplest case of 
photoluminescence; it is called "resonance radiation". However, if, as indicated in Figure 
2, several other levels C, D, ... are located between N and F, other transitions to C, D, ... 
can occur, resulting in the emission of spectral lines of frequency Vrc, Vrp, ... These 
frequencies must be smaller than Vrw and therefore photoluminescence can have no 
greater frequency or shorter wave-length than the exciting light [Ref. 11: p. 2] 

1. Perturbation Theory 

For a molecule that absorbs or emits quanta of radiation under the influence of an 


electric field, perturbation theory of quantum mechanics provides the theoretical model 
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1. Resonance Radiation 


2. Phosphorescence 


3. Fluorescence 
4. and 5. Anti-Stokes Fluorescence 

Figure 2. Energy Level Diagram for the Representation of Fluorescence and 
Phosphorescence? 





that defines the physical interaction. The general solution ofthe wave equation for a 
molecule in the presence of a perturbing field results ın certain derived relations that 
govern the changes of the molecule. One such outcome mandates that for radiation to be 
absorbed and transition from a stable ground state of energy to an electronically excited 
state of energy to occur, the frequency of the quantum must be proportional to the 
difference in energy between the states, a situation previously discussed. Once this 
condition is met, the probability of the transition is governed by the value of the transition 


moment integral, a second corollary from the solution of the perturbed wave equation: 


“Figure 2 is found in Ref. 11, page 4, as Figure 2 


т = Y, MY, dí 

Неге, Ғ, апа“) аге (һе total wave functions for the upper and lower states 
respectively ( the wave functions result from the unperturbed solution of the Schrodinger 
wave equation ) and M 1s defined as the dipole moment operator [Ref. 10: p. 583]. 
Evaluation of this integral determines the probability of an optical transition between two 
states. Detailed T— has shown that transitions are probable only between states 
characterized by wave functions of a certain class. The results have been summarized in 
the form of selection rules which are formulated in terms of the quantum numbers that 
describe the electronic states involved [Ref.12: p. 51]. 

Perhaps the most important is the rule governing spin multiplicity. This states that 
the spin of an electron must not change during an electronic transition. As will be seen 
later, although from a theoretical standpoint transitions between states of different 
"multiplicity" do not intercombine, in practice they are observed, but the probability is 
usually much lower than an allowed transition. Thus, selection rules are good empirical 
guides, but they are not inviolable.[Ref. 13: p. 36] 

Although the selection rules considered have been based on the absorption and 
emission of light by a molecule under the influence of a perturbing field, a second set of 
selection rules similar to the first but derived under a different set of circumstances govern 
radiationless processes. These processes, in addition to the spontaneous emission of light, 
compete with the excited state of a molecule in depleting its acquired energy and allowing 


it to return to its ground state. 


2. Spontaneous Emission 

An excited state may undergo a third radiation-related transition to a lower state by 
means of spontaneous emission, 1.e., fluorescence and phosphorescence. In this process, 
an excited species looses energy in a spontaneous manner in the absence of an electric 
field. Itis a random process: the rate-of-loss of excited species following a kinetic 
first-order reaction. Ifa process follows first-order kinetics described by a rate constant 
"К", the mean radiative lifetime, (the time t taken for the intensity to fall to 1/e of its initial 
value), is given by 

t= 1/k 
When a radiative transition occurs between states of similar multiplicity, it corresponds, 
according to the selection rules, to an allowed transition. Such transitions are 
characteristic of fluorescence emission having lifetimes equal to t = 10°® seconds. 
However, there are transitions in which t is much higher, 10° seconds. When selection 
rules suggest a low probability for a transition to occur, such as between states of different 
multiplicity, the transition, though theoretically forbidden, can occur, but the process 1s 
relatively slow. This 1s representative of phosphorescent emission. [Ref. 10: p. 585] 
3. Triplet and Singlet States 

In the absence of a magnetic field, the energy of a single unpaired electron 
characterized by its given values of "n" the principle quantum number) and "I" (the 
azimuthal quantum number) must be further refined to explain the observed fine structure 


of atomic spectra. The additional amount of energy that the spectroscopic evidence 
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indicates is attributed to the interaction that occurs between the angular momenta from the 
intrinsic spin and orbital motion of the electron. When several electrons within an atom 
(or molecule) are present, the resultant orbit and resultant spin angular momenta 
components can combine in several ways resulting in several energy levels for a single 
configuration of an atom (or molecule). This interaction is known as Russell-Sanders 
coupling, the result of which is the "multiplicity" of energy levels that a given electron 
configuration of an atom may possess. Thus, a single energy level for a given electron 
configuration for an atom (or molecule) is referred to as singlet state, two energy levels as 
a doublet, three energy levels a triplet, etc. [Ref.14: p. 79] A more thorough discussion is 
given in Appendix A. 

Nearly all molecules have a ground-state electronic configuration in which the spin 
of each electron is opposed to, or paired up with, the spin of another electron which 
otherwise has the same spatial quantum numbers, a consequence of the Pauli exclusion 
principle. In this instance the total angular momentum of each electron within the pair, 
composed of the spin and orbital momenta, tend to cancel with each other . When all the 
pairs of electrons within the molecule are considered, the net result is zero total angular 
momentum (vector), producing a single energy level for the given electron configuration. 
Molecules will have, in addition to a ground state that is usually a singlet state, a number 
of excited states that also have all the electrons paired and that are, therefore, also singlet 


States. [Ref. 15: p. 16] 
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Whether or not the ground state 1s a singlet state for molecules with an even 
number of electrons, there will be excited states in which a pair of electrons have their 
spins in the same direction, giving the molecule a net total angular momentum (vector). 

In this case, the coupling interaction results in three energy levels for the configuration, 
and thus the designation triplet state. 

Therefore, promotion of an electron to an upper state can in principle give rise to 
two types of energy states, the triplet and singlet, while still maintaining the same electron 
configuration. Therefore, for every singlet state there is a corresponding triplet state, with 
the triplet state displaying a lower energy level. Figure 3 illustrates the two types of 
electronically excited states, where the spin orientations for each state are depicted as 
arrows next to the appropriate state. The additional lines within each electronic state refer 
to the various vibrational and rotational energy levels the molecule may acquire with the 
absorption of energy. 

4. Intramolecular Processes Involving Excited States 

With the absorption of a quantum, a molecule is raised from a singlet ground to a 
singlet excited state. The selection rules stipulate that the absorption of electromagnetic 
radiation does not "unpair" the electrons of a molecule (transitions =н between states of 
like multiplicity) thus preventing direct absorption of light from a singlet to a triplet state. 
The molecule, having been raised to an upper vibrational level within the excited state, 


begins to lose the excess vibrational energy by collisions with surrounding molecules. The 
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Figure 3. Energy Levels of Molecular Excited States and Transitions Between Them? 
radiationless process of internal conservation then occurs whereby the molecule passes 
from a low vibrational level of the upper state to a high vibrational level of a lower state 
having the same total energy. Once internal conversion is complete, this cycle of events 
repeats itself leading to the net result of the molecule falling to the lowest vibrational level 
of the first excited singlet state as illustrated in Figure 3. [Ref. 15: p. 9] 


At this point, the molecule may return to the ground state by one of three paths: 


S1 > So (Internal conversion) 
Sı > So + hv (Fluorescence) 


S; — T; (Intersystem crossing) 


? Figure 3 is found in Ref. 16, page 71, as Figure 5-1. 
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The internal-conversion process, involved with the intramolecular energy transfer between 
the first excited singlet and ground state, takes place at a much slower rate (10? seconds) 
as compared with the same process occurring between the upper excited singlet states to 
the first excited singlet state (10 seconds). The increase in inefficiency in energy transfer 
is explained by the relatively large energy separation between S, and S, as compared to 
the energy differences between the upper excited states. This allows fluorescence (10° to 
10" seconds) and, with the combined effect of the spin-forbidden nature of the states, the 
process of intersystem crossing (10? seconds) to effectively compete with the internal 
conversion process. 

When a molecule returns to the ground state via the first excited triplet state, it 
undergoes a conversion process known as intersystem crossing. Intersystem crossing is a 
non-radiative process between states of different multiplicity. It is theoretically forbidden 
with regard to the formal selection rules, but in practice, these transitions do take place 
(although with extremely low probability). This is made possible through the 
"breakdown" of the spin-orbit coupling model discussed in Appendix A. In the absence 
of magnetic or electric fields the total angular momentum of the molecule remains 
constant, and the ability to discern between specific states such as singlet and triplet is 
possible. However, when an electromagnetic field is present, though the total angular 
momentum vector remains constant, the spin and oribital degrees of freedom become time 
dependent, and it becomes meaningless to speak of specific states [Ref. 16: p. 12]. This 


ultimately leads to the mixing of states. The interaction is enhanced by heavier elements 
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due to an increase of the magnetic field developed by the nucleus, as in the case of 
platinum octaethlyporphyrin (PtOEP). 

The lowest vibrational level of the lowest triplet state is generally situated at an 
energy level below that of the lowest excited singlet state, but its upper vibrational levels 
reach to the bottom of the singlet level, and intersystem crossing occurs with the molecule 
crossing over to one of these upper vibrational levels of the lowest triplet state. From here 
the molecule rapidly looses its excess vibrational energy and falls to near the lowest 
vibrational level of the triplet state. 

Again, the molecule may return to the ground state from this point by one of the 


following paths: 


Tı —S. (Intersystem Crossing) 


T; — S, * hv (Phosphorescence) 


As before, it might be expected that the analogous process of intersystem crossing from 
the lowest vibrational level of the triplet state to the ground state would take place with 
equal facility as the transition from S; — T,;. In fact it is usually some 10° - 10° times 
slower. If radiationless T, — S,transitions occurred with rate constants of 10* seconds, 
the molecule would have no time to emitT; — So phosphorescence by the spin-forbidden 
radiative process. On the other hand, if the S1 — T, intersystem crossing process were as 
slow as that of the T; — S, radiationless process, then it could not compete effectively 
with fluorescence. Phosphorescence would then be a rare phenemenon and the 


photochemical behavior would be very different. The difference between the rates of the 
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Sı — Tı and Tı > Soradiationless processes is, as was the case with internal conversion, 


a function of the energy differences between S, and Т, , and, T, and S,.[Ref. 15: p. 41] 


D. OXYGEN QUENCHING 


The exceptionally long radiative lifetimes of molecules in the lowest triplet state make 
them vulnerable to a variety of radiationless processes that deactivate the molecule, 
allowing it to return to the ground state. As an example, they are very susceptible to 
collisions with molecules, and in such cases may encounter molecules of a specific nature 
which are particularly effective in removing excitation energy. 

Previous reference was made to the fact that most molecules with an even number of 
electrons have a ground state in which their electrons occupy orbitals in pairs. These 
molecules were said to have singlet states. There are however some molecules containing 
an even number of electrons for which the states of lowest energy have two unpaired 
electrons occupying different orbitals. The most important of these is molecular oxygen, 
of which the ground state is a triplet. Molecular oxygen is very effective in removing 
triplet excitation energy, and can act at exceedingly low concentrations. [Ref. 15: p. 36] 

The energy transfer that takes place when an excited triplet PCOEP molecule collides 
with a triplet oxygen results in a ground-state singlet PCOEP molecule and a singlet 
excited oxygen molecule. The reaction is expressed as follows: 

Ti *?0O; 9S, *!O; 
The reaction occurs very rapidly but is consistent with the formal selection rules. It is this 


bimolecular reaction that results in the quenching of the phosphorescence. The University 
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of Washington found that when PtOEP was embedded in an oxygen-permeable silicon 
resin and subject to atmospheric pressure, the quenching reaction had a rate comparable to 
that of phosphorescence. Variations in atmospheric pressure from 0 to 1 atmospheres 
demonstrated a suitable drop in phosphorescence yield of PtOEP that would allow it to be 
useful as a measurement tool in the study of static pressure changes in aerodynamic wind 


tunnel testing. [Ref. 9: p. 18] 


E. MATHEMATICAL MODEL DEVELOPMENT 


Oxygen quenching of luminescence by collisional deactivation of the emission process 
is the basis for developing a prediction model. It is important to note that the end result of 
this development, the Stern-Volmer relation, is strictly valid only for free molecules in 
solution. However, it has been found that it describes accurately the quenching process 
that occurs in the present situation. The basis for the present development is the work of 
Parker [Ref. 15], which may be consulted for further details. 

The development first begins with the definition of the triplet formation efficiency, $. , 
being the number of triplet molecules formed per quantum of exciting light absorbed. 
Assuming now that a system, illuminated with a beam of exciting light of constant 
intensity, is in a photochemically steady state. The rate of population of the lowest excited 
singlet state, S, , 15 found to equal the rate of light absorption, I, , and the rate of formation 
of triplet molecules by crossing from S, is equal to I. . 

Figure 4 illustrates the rates of the various processes by which the triplet is produced 


and consumed; these are represented in equation form as follows: 
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Figure 4. Summary of Transition Types and Associated Rates* 


S. 381 Т! Rate = 1,6, (Triplet Production) 
Tı — So + ћу k [T] (Phosphorescence) 

Т, > S, + heat km[T] (Intersystem Crossing) 
Tit+tq->Sotq kg[T][q] (Triplet Quenching) 
Т, шаш > $; k.[T] (Intersystem Crossing) 


The square bracketed [q] and [T] represent concentrations. [q] represents the 
concentration of all quenching impurities in the system, k,{q] a composite pseudo 
first-order rate constant for impurity quenching, and k.[T] represents the rate of thermal 


activation of triplet molecules back to the excited singlet level [Ref. 15: p. 87]. 


*Figure 4 is found in Ref. 15, page 69, as Figure 23 
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[T] represents the concentration of excited triplet molecules when the steady state has 
been achieved. All processes so far considered in Figure 4, with the exception of light 
absorption, are first-order reactions having rate constants equal to the reciprocals of the 
lifetimes quoted. 

Now consider, in addition to the previous statements and definitions made with respect 
to the steady state of the system, that the rate of production of excited triplet molecules 
(in stateT, ) is just balanced by their rate of disappearance via the various processes 
discussed above. This can be expressed as 

Гаф: = (kp + km + kgfq] +ke){T] (1) 
Since the rate of emission of phosphorescence is k,[T] , the phosphorescence efficiency 
Op is given by 
p = kp[T] / Ia (2) 
Equation (1) may be solved for I, and substituted into equation (2) to arrive at 
фр 7 Крфі / (Кр + Km + kg[q] + ke) (3) 
In the absence of quenching agents affecting the removal of triplet excitation energy, 
equation (3) gives, for the phosphorescence efficiency, 
фр = Коф. / (Кр + Къ + Ко) (4) 
Taking the ratio of equations (3) and (4): 
фр/фр = 1 +kg[q] / (Kp + Km + Ko) | (5) 


Ог, 
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Фр/фр — 1 * КІ] (6) 
where 
K 7 t,k, and 1/т, = (К + К. + К.) (7) 
Equation (6) is known as the Stern-Volmer equation for quenching, and the constant K as 
defined in equation (7) is referred to as the Stern-Volmer quenching constant. 
For this particular application, the quenching impunity, [q], may be defined as 
[q] =S Po, (8) 
where S is the solubility of oxygen in the system and Po, is the partial pressure of oxygen 
in equilibrium with it. This leads to the relation 
ф/ф%,-1%К Ро, (9) 
Finally, by definition, the luminescence, 1, is linearly proportional to the quantum yield 
( 9 = photons emitted/photons absorbed) for a given quenching condition, (i.e.,] = Las ) , 
where L,. 1s the excitation light intensity by the luminescent species. Thus 
po (10) 
where Iaa is the maximum possible value that occurs in the absence of a quencher. 
1. Temperature Effects 
The effect that temperature has on the efficiency of phosphorescence may lead to a 
decrease in the intensity by several orders of magnitude over the range of 77° K to room 
temperature. The effect results in the long lifetime of triplet states and the 
correspondingly greater susceptibility to impurity quenching and intersystem crossing 


[Ref. 15: p. 89]. 
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The rate constants for each of the intramolecular processes depend on temperature, 
following the Arrhenius equation, 
Кае 
where k is the rate constant, A is the pre-exponential factor, E, is the activation energy, R 
is the gas constant, and T is the temperature in degrees Kelvin [Ref. 9: p. 20]. Hence, in 
equation (10), I, I_,. , and K are all functions of temperature. 
2. Aerodynamic Application 

In order to apply the form of the Stern-Volmer relation presented in equation (10), 
it would be necessary to purge the vicinity of the surface of interest of all oxygen in order 
to determine І, , . In most cases of experimental work, this is impractical. A more 
convenient approach is to manipulate the governing equation so that determination of I, ,. 
becomes unnecessary. 

Equation (10) is first modified by noting that the oxygen partial pressure, Po, , is 
related to the air pressure, P, by Po, ^ X P, where X is equal to the mole fraction of 
oxygen. Equation (10) then becomes: 

Imax/I=1+K’ XP (11) 

In the method used in the present work, two images are taken. In the "flow-on" 
case, the pressures are unknown. In the "flow-off" case, the pressure distribution is 
constant on the surface and equal to atmospheric pressure. These conditions are 


represented, respectively, as follows: 
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Г.и ТЫ ӘР (12) 
and ET, = Kae (13) 
where the subscript "o" denotes the value for the flow-off conditions. Dividing equation 
(13) by equation (12), 
1Л=(1+К ХР)(1+К' ХР,) (14) 
Rearranging, 
1,Л-1/41-К ХР,)-((К ХР,У(1-К ХРОМРР,). (15) 
which gives 
I/A = A +В (РЉ,) (16) 
where, 
А=1/(1+К ХР„)ап4В=(К ХР„)/(1+К ХР.) (17) 
A and B represent the coating sensitivity coefficients and are, in general, functions of 
temperature for both wind-on and wind-off conditions. In the absence of a temperature 
variation across the surface, A and B have the same value over the whole image. Notice 
that A+ B=1. This condition occurs only when the images of I and I, are taken at the 
same temperature. If any temperature variation is present between the two images, a 
value other than one will result. [Ref.1: p. 34] Techniques discussed in Section 3 will 


outline the procedures used for determining these coefficients. 
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III. EARLY DEVELOPMENT USING A FREE JET APPARATUS 


A. APPARATUS AND FLAT PLATE ASSEMBLY 
A sonic-jet nozzle system located in the Gas Dynamics Laboratory (Bldg. 216) at the 
Naval Postgraduate School was used in the initial experiments. A schematic of the facility 


is illustrated in Figure 5. A photograph of the jet nozzle is shown in Figure 6. 
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Figure 5. Schematic of the Underexpanded, Sonic Jet Nozzle and Supporting 
Facility 
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Figure 6. Underexpanded Jet Nozzle Apparatus 


The flat plate used in the present work was a model designed and previously used for jet 
interaction studies in the supersonic wind tunnel. The plate was made of stainless steel 
with a leading edge diameter of 0.008 inches. Twelve of the existing pressure taps were 
chosen as close to the centerline of the plate as possible. Figure 7 shows a schematic of 
the surface of the plate showing the twelve taps relative to the centerline. The 
pre-arrangement of the taps, in addition to discovering two plugged ports, required that 


three offset taps be used to maintain near uniform spacing in the streamwise direction. 
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Figure 7. Schematic of the Pressure Port Arrangement on the Plate Surface 

A simple bracket was made to mount the model so that the surface of the flat plate 
become a radial plane, effectively dividing the ax1-symmetric jet in half. The bracket and 
model are shown in Figure 8 while Figure 9 illustrates this assembly mounted to the face 
of the jet. Plate washers were placed on the mounting bolts to create a 1/4 inch gap 
between the leading edge of the plate and nozzle face. This was necessary as it was found 
that the jet would not properly form if the plate were installed with its leading edge at the 
plane of the nozzle opening. Pneumatic pressure measurements were made using the 


Scanivalve ZOC system developed by Wendland [Ref. 17]. 
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Figure 9. Bracket and Model Assembly Mounted to the Face of the Sonic Jet 
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1. Freestream Flow Field of the Underexpanded Sonic Jet 
The mechanics of formation and structure of an underexpanded jet from a sonic 
nozzle have been studied in detail by Adamson and Nicholls [Ref. 18: p. 265]. A 


schematic of the structure is illustrated in Figure 10. 
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Figure 10. Sketch of Jet Structure Behind a Highly Underexpanded Nozzle” 
The "Mach disc" within the free-jet offered an event that would produce a pressure 
discontinuity similar to that of the starting-shock structure in the supersonic wind tunnel, 
but without the visual restrictions imposed by thick test-section windows. 
Using the schlieren technique, visualization studies of the flow field across the flat 


plate were recorded. The data acquired using the pressure-sensitive paint on the flat plate 


"Figure 10 is found in Ref. 18, page 265, as Figure 1 
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could then be interpreted with reference to the flow structure. These studies began by 
observing the underexpanded jet without the flat plate mounted to the jet nozzle. This 
was done to provide a reference for changes in the flowfield as a result of the interaction 
with the plate. Since the jet structure was axi-symmetric, it was anticipated that, with the 
plate surface aligned along the centerline of the jet, the plate would simply isolate half of 
the jet structure. Figure 11 illustrates the jet ана without the flat plate installed at 65 


psig plenum pressure. The flow proceeds from right to left. 





Figure 11. Underexpanded Free Jet at 65 psig Plenum Pressure 


The flat plate was then bolted to the nozzle face and several runs were made at 
varying plenum pressures to record the shock structure. Three specific plenum pressures, 
65, 80, and 90 psig were chosen as test points at which to compare pressure 
measurements with data using the luminescent paint, and the schlieren images of the flow 
field. The schlieren images of the three events and the corresponding pressure 


measurements from the orfices along the flat plate are illustrated in Figures 12, 13, and 14. 
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Figure 12. Flow Visualization and Pressure Distribution Across Flat Plate at 65 psig 
(Flow is from left to right) 
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Figure 13. Flow Visualization and Pressure Distribution Across Flat Plate at 80 psıg 
(Flow 1s from left to right) 
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Figure 14. Flow Visualization and Pressure Distribution Across Flat Plate at 90 psig 
(Flow is from left to right) 
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The images and the pressure measurement graphs have been aligned so that the pressure 
changes seen on the graphs correspond to the events occurring on the plate surface. 

As can be discerned from the schlieren images in Figures 12, 13, and 14, two 
oblique shocks appear, a consequence of the plate being placed in the flow field. Figure 
15 shows a comparison of the flow fields with and without the plate. The oblique shock 
emanating from the leading edge of the plate can be explained as being due to the 


underexpanded process the flow undergoes as it leaves the nozzle. The flow, as it exits 





Figure 15. Flow Field Comparison Without (Top) and With (Bottom) the Flat Plate 
the nozzle, is off-axis with respect to the centerline of the plate thus encountering the plate 


at an incidence angle other than zero degrees. The flow in these images 15 from right to 


left. 
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The second oblique shock's structure gives the perception that 1t extends from the 
surface of the plate to the jet boundary; however, closer examination of the two images in 
Figure 15 shows that the oblique shock actually merges with the reflected shock (see 
Figure 10) giving the erroneous illusion of a single entity. Of greater interest is the extent 
of the oblique shock to the point of total suppression of the normal shock. The sensitivity 
of the flow to the plate was much greater than first thought, and 1t was unclear what 
events were taking place to justify the magnitude of the oblique shock. Additional study 
seemed warranted. 

a. Further Flow Field Studies 
To visualize the flow, methanol was injected through the second pressure port 
in the plate surface. Several light fixtures were placed around the plate to better illuminate 
the injection area and to provide contrast between the background and flow field. A video 
recorder was used to record the events. 

Figures 16 through 18 show the sequence of events from selected frames of 
the video at two different vantage points. In each set of images, the flow proceeds from 
left to right. 

In Figure 16, the plenum pressure was steady at 30 psig -- the boundary layer 
up to this point remained attached. At 38 psig, Figure 17 shows initial indications of 
reversed flow at the center of the plate. In addition, 1/2 inch to the left of this point, a 


pooling of methanol formed. This behavior indicated the presence of a bubble in which a 
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Figure 16a. Flow Visualization at 30psig -- Cross Flow View 





Figure 16b. Flow Visualization at 30 psig -- Acute View 
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Figure 17b. Flow Visualization at 38 psig -- Acute View 
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Figure 18a. Flow Visualization at 65 psig -- Cross Flow View 





Figure 18b. Flow Visualization at 65 psig -- Acute View 
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portion of the injected mass would be expected to stagnate near the separation point, as 
the flow re-circulated within the bubble's boundaries. This is more clearly shown in the 
schematic of Figure 19. With increasing plenum pressure, the separation and recirculation 
region of the flow field gradually moved downstream, eventually encompassing the 
injection port, and allowing the alcohol to more clearly define the boundary of what was 


indicated to be a separation bubble. 
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Figure 19. Details of Separation Bubble 
To verify whether flow reattachment was in fact occurring downstream, a 
second injection port, port number 5, was rigged to concurrently admit methanol into the 
flow. Figures 20 through 23 show the sequence of images recorded in one such test. 
Figure 20 shows that at 25 psig, attached flow was prevalent. At 40 psig, in 


Figure 21, the injectant at the first port (port number 2) indicated the presence of reversed 
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flow, while the injectant at the second port (port number 5) continued to indicate attached 
flow. The vortex patterns are an indication of the pressure gradients that exist within the 
bubble. In Figure 22, with the plenum pressure at 68 psig, the first sign of reversed flow is 
evident at the second injection port, while Figure 23 shows the port fully immersed in 
reversed flow at 75 psig. 

The onentation of the plate relative to the centerline of the nozzle exit, and 
thus the centerline of the flow, was found to be 89.7 degrees. Since alignment could 
significantly influence the flow characteristics, the plate's inclination was changed to 90.3 
degrees and the test runs were repeated. No dissimilarities between the results obtained at 
corresponding plenum pressures were noted. 

The evidence of reversed flow and reattachment strongly suggested the 
presence of a closed separation bubble within the underexpanded jet when the plate was 
present. The apparent extent of the bubble, reaching 3/8 to 1/2 inch above the plate, is 
very surprising. However, it is consistent with the oblique shock in place of the normal 
shock. Further data are needed, such as impact probe measurements, in order to fully 
define the flow structure. For the present purpose, which was the investigation of 


pressure sensitive paint, this was not needed. 


38 





Figure 20. Flow Visualization at 25 psig 





Figure 21. Flow Visualization at 40 psig 
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Figure 23. Flow Visualization at 75 psig 
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B. OPTICAL MEASUREMENT SYSTEM 

The measurement system was composed of three main components: (1) the 
luminescence coating, (2) the illumination source, and (3) the luminescence detector and 
associated data processing hardware and software. Each component is discussed 
individually below. A schematic of the system used in conjunction with the free jet 
apparatus is illustrated in Figure 24. 

1. Luminescence Coating 

The luminescence coating used in the present work was developed by the Chemistry 

Department at the University of Washington and was provided, as a sample, by the 


NASA-Ames Research Center. 
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Figure 24. Schematic of Pressure Measurement System 
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The coating contained platinum octaethylporphyrin (PtOEP) as its active molecule. 
PtOEP has characteristics of high luminescent quantum yield (~90%), short triplet life 
(~100 micro seconds), and a high sensitivity to oxygen quenching 
[Ref. 1: p.35]. Figure 25 shows the excitation and emission spectra for the molecule. 
Excitation (absorption) peaks at 380 or 530 nm, whereas luminescence peaks at 650 nm. 

The coating mixture consisted of the active molecule dispersed in a Genessee 
GP-197 polymer-resin solution. The concentration was 0.1 gm of PtOEP per liter of 
polymer resin [Ref. 19]. When the mixture was applied to the surface, the solvent 
evaporated leaving a hard, thin film of PtOEP dissolved in an oxygen-permeable polymer 
[Ref: 9: p. 35]. Several undesirable characteristics associated with this coating are 


described in Appendix B. 
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Figure 25. Excitation and Emission of Platinum Octaethylporphyrin 


42 


a. Reflective Backing 

Before applying the luminescent coating to the surface, a white backing was 
first applied. A significant improvement in luminescence performance was thereby gained 
[Ref. 9: p. 32], allowing the camera to see a larger signal and an enhanced signal to noise 
ratio. 

The particular type of white paint used for this backing was an important 
consideration. Kavandi [Ref. 9: p. 33] found that some paints tend to darken after long 
exposure to air and ultraviolet light. Krylon glossy-white spray paint (#1501) reportedly 
offered the greatest resistance to this effect and was therefore used in the current work. 

2. Illumination Source 

McDonnell-Douglas [Ref. 6: p. 3] considered several types of illumination sources 
to provide an adequate number of photons within the appropriate excitation wavelength of 
the paint. The results showed that tungsten-halogen incandescent sources were the most 
suitable in terms of compatibility and performance. Although the data in Figure 25 would 
suggest that arc lamps emit more energy in the excitation range of the coating, they also 
emit an EMF pulse that is damaging to CCD arrays. In addition, tungsten-halogen lamps 
are extremely stable, with minimum intensity fluctuations, and any such variations will 
directly effect the accuracy of the results. Consequently, an Oriel Corporation 1000 W 
quartz tungsten-halogen lamp, Model 66200, and the associated controller, Model 6405, 


were selected for the present measurement system. 
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Spatially-uniform illumination was also essential to reduce measurement errors. 
(According to McDonnell-Douglas [Ref. 6: p. 4], this 1s especially true when the Е 
are highly curved). Non-uniformity was found when the published procedures 
[Ref. 20: p.13] were followed to produce a collimated beam from the light source. The 
| procedure resulted in generating a focused image of the filament in the lamp. However, 
by reversing the orientation of the barrel containing the lenses, a greater range of 
adjustment was achieved. Figure 26 illustrates this modification, which made it possible to 
produce a more uniform beam, although a slight degradation in spectral intensity was 


noted. 
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Figure 26. Modification of Focusing Assembly of the Light Source 
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a. Filters 

The use of a broad-band light source mandated the use of filters to ensure 
that the appropriate wavelength for excitation of the luminescence was allowed to pass 
while blocking unwanted light. Unwanted light was that which would pass through the 
narrow-band filter that was installed on the camera. 

As Figure 25 shows, there are two wavelengths at which excitation of the 
luminescent paint is optimum. Although either wavelength could have been used, the 
380 nm wavelength was chosen because of its greater separation from the emission 
wavelength of 650 nm. This also had the advantage, as discussed in [Ref. 6: p. 5], of 
allowing a filter with a more gradual cutoff to be used which would allow more light 
energy to be utilized to excite the paint. 

A combination of filters was used in the illumination system to isolate the 
380 nm wavelength. The source beam first passed through a dichoric or cold mirror, Oriel 
Model 66228, which had a spectral range of 350 - 450 nm and a cut-off frequency of 
550 nm [Cold mirrors are oriented at 45 degrees relative to the incoming beam from a 
source. The infrared wavelengths are transmitted undeviated, while the ultraviolet 
wavelengths are reflected through 90 degrees. This prevents the infrared from being 
reflected back to the source and overheating the system. ] 

The source light was then directed through an Oriel interference filter, Model 


$7521, having a bandwidth of 70 nm centered at 400 nm. The filters were installed in a 90 
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degree beam turner and mirror holder, Oriel Model 66246, which was fitted to the lamp 
housing. 
Finally, the camera used to record the emitted luminescent energy from the 
paint had a 650 nm interference filter with a 70 nm bandwidth, Onel Model 57610, 
mounted over its lens to block out all other frequencies. 
3. Detection System 

Luminescent data were acquired using a COHU Model 4910 high performance, 
monochrome, low-intensity CCD camera and transferred to a framegrabber board in a 
80386 PC using an R-59 coaxial cable. The camera was fitted with a Cosmicar 16 mm, 
f/1.4 lens. To ensure the output signal would be linearly proportional to light intensity 
received, the camera was configured with the gamma set equal to one and the automatic 
gain control (AGC) disabled [Ref. 2: p. 3342]. The camera's spectral response, as shown 
in Figure 27, shows that its relative response at 650 nm is approximately 65%. 

The camera produced a frame every 1/30 of a second and thus was limited to 
recording pressure fluctuations at rates less than 30 Hz. [A simple test to determine if the 
illumination system was bleeding over into the detection system was to illuminate an 
unpainted surface (1.e., no luminescent paint) and then direct the camera towards this 
surface to see if any spurious light was being detected. None should be present]. 

The analog output signal from the camera was initially digitized using an IDEC 
Supervision 16 framegrabber board residing in a 80386-based PC. The raw data were 


then reduced and analyzed using "Image Pro Plus" analysis and processing software. 
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Figure 27. Spectral Response Curve for the COHU cM Model 4910 

However, this hardware/software combination was extremely limited in application 
resulting in long data acquisition periods and time-consuming data reduction. Because of 
the characteristic degradation of the luminescent coating, the accuracy of the data was 
questionable. 

The initial hardware/software combination was replaced with a much more versatile 
Epix 4MEG Video Model 12 framegrabber board and software which included image 
analysis and processing programs. The spatial resolution of the framegrabber was 
operator adjustable, thus a balance between sampling resolution and the board's memory 
storage could be tailored to a given application. For the underexpanded-jet and fence 


work, the spatial resolution was set to 752 x 480, the highest resolution setting. The 


gray-level resolution of the board was 8-bits. 
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C. EXPERIMENTAL PROGRAM 
1. Flat Plate Preparation 

The surface of the flat plate was cleaned with acetone, ensuring the removal of all 
oils and dirt. Fine wire was then inserted into each of the pressure ports to prevent paint 
from entering the holes. The white Krylon undercoating was then applied in a series of 
three to four thin coats, allowing four to five minutes between coats. Ten minutes after 
applying the last coat, the fine wires were removed from the ports. This was necessary to 
prevent the paint from drying on the wire. [Significant portions of the undercoating could 
be removed from around the pressure port when extracting a wire to which the paint had 
hardened]. At this point, the white undercoating was allowed to dry for one hour. 

Prior to the application of the pressure-sensitive paint, the fine wires were 
re-inserted into the ports. The luminescent coating was then applied to the surface using a 
standard hobby airbrush. The nitrogen supply in the Gas Dynamics Laboratory was used 
as the pressure source for the airbrush (instead of shop air) as it offered a contaminant-free 
pressure source. With a properly adjusted airbrush and a good painting technique, 
continuous application of the luminescent coating to the surface could be achieved until 
the surface was "suitably" covered. Although Kavandi et al [Ref. 2:p. 3343] refers to 
normal coating thickness as being in the range of 5 to 15 micrometers, it was both 
impractical and unnecessary to measure the coatings to determine if the paint fell within 
these prescribed limits. The approach used in the present work was simply to observe 


when the pigment of the pressure-sensitive paint seemed insensitive to further application, 
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and that the texture of the surface appeared smooth. Several attempts were necessary to 
develop the technique and the judgement in applying the luminescent coating. Three to 
four minutes after the completion of the coating, the fine wire were removed. Once the 
model was installed for testing, sufficient time had elapsed for the coating to have properly 
cured. Note that it was important to avoid touching the painted surface at any time since 
body oils could adversely affect the behavior of the luminescent coating. 
2. Experimental Procedure 

Setting the video resolution for the framegrabber to the maximum (752 x 480) 
resulted in 11 frame buffers being available to store for image data. The controlling 
software for the framegrabber was capable of capturing a sequence of images, digitizing 
them, and storing each image of the sequence in consecutive image buffers, the number of 
buffers used in the sequence (from one up to the number of available buffers) being 
specified by the operator. Because of the limitation of collecting only 11 images for each 
execution of the sequence algorithm, and the fact that averaging one hundred images was 
desirable to obtain noise-free data, for each of the wind-on and wind-off conditions, ten 
runs were executed. (10 frame buffers vice 11 buffers were utilized). Since the 
underexpanded free jet was a steady and repeatable event, this method of collecting data 
was a feasible solution for the limited memory of the frame-grabber board. 

The data for the wind-off condition was collected first. The flat plate was 
illuminated with the tungsten-halogen light source while the laboratory lights were turned 


off, preventing any spurious wavelengths from entering the camera and offsetting the 
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image data. A real-time image of the plate was displayed on the image monitor so that the 
camera Output could be adjusted for the best quality picture using a combination of focus, 
manual gain, and lens iris. Manual gain was necessary because of the low light condition, 
but care had to be exercised because with an increase in gain there was an accompanying 
increase of noise introduced into the image. The shutter speed for the camera was set to 
off. 

After the first sequence of images was captured and stored in the frame buffers, the 
light source was switched off to minimize the degradation to the paint. The light source 
was such that the line voltage applied to the lamp housing was adjustable; however, power 
was ultimately controlled with an on/off switch. Thus, the same illumination intensity was 
maintained each time the lamp was turned on. 

The ten images were averaged together to form a single image. The single image 
was saved to a floppy disk. The light source was again switched on in preparation for the 
second sequence of images, allowing time for the luminescent coating to reach its 
maximum intensity. (This is referred to as the "induction period", which is explained in 
Appendix B). The second sequence of images was then captured and stored, and the 
tungsten-halogen lamp shut-off. The images were averaged to give a single image and 
again saved to floppy disk. This routine was executed until ten averaged images for the 
wind-off condition had been saved to floppy disk. Each sequence of capturing the data 


took less than one minute. 
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The procedure for obtaining images for the wind-on condition was similar. With 
the light source on and the plate properly illuminated, the underexpanded-jet apparatus 
was started. Once the desired pressure was set using the plenum control valve, and the 
pressure was stabilized, the imaging sequence was executed. When the software had 
completed its routine, the jet was shutdown and the light source turned off. The jet 
apparatus was equipped with a shutoff valve so that the jet could simply be turned on and 
off without disturbing the setting of the plenum control valve. This ensured repeatability 
of the event for each sequence of 10 images. 

It is important to note that throughout these events, the camera and light source had 
to remain fixed, spatially, with respect to the plate. In addition, no adjustments of any 
kind could be made to the instruments or equipment. Any deviations would introduce 
errors into the results. 

Conventional pressure data were obtained concurrently with the luminescent paint 
data using the CALS YS2000 data acquisition system [Ref. 17]. Twenty samples per port 
using a data collection rate of 33 Hz were obtained and averaged together. These pressure 
data would be used later in the calibration process. Ambient pressure and temperature 
were also recorded. 

3. Image Data Reduction 
a. Image Processing 
Several image-processing techniques were used to manipulate the data in 


order to obtain the desired result. The techniques are described individually below. 
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Appendix C provides a summary of the steps used in operating the EPIX software for 
acquiring, processing, analyzing, and reducing the data. 

Round-off error from the image processing routines was a concern. To 
ensure against loss of accuracy due to round-off error, the system hardware had to be 
capable of doing at least 16 bit arithmetic. The EPIX frame grabber board was capable of 
conducting 32-bit operations. When the image was processed, the 8-bit intensity values 
for each pixel in the image were first converted to 32 bit words. Once the processing was 
complete, the modified values of each pixel were converted back to an 8-bit format for 
display. If the mathematical operations were conducted with anything less than 16-bit 
precision, all significant digits would be truncated. 

(1) Averaging. The emission of photons from any source is a random 
process following a Poisson distribution, and the number of photogenerated carries 
collected in a potential well of the camera array (see Appendix D for an explanation on the 
operation of a CCD camera) in time "t" is thus a random variable. The standard deviation 
of this random variable is referred to as photon (shot) noise. This noise is a fundamental 
limitation in all imaging applications since it is a property of light itself, not of the image 
sensor. 

By averaging a sequence of images, the signal-to-noise ratio of the data 
has been found to be significantly improved. This process has the effect of reducing the 


scatter in the luminous intensity measurement [Ref. 6:p. 8]. McDonnell-Douglas [Ref. 6] 
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has found that as few as eight images averaged together showed a marked reduction of 
noise on the signal. 

The reason for choosing one hundred images for each of the wind-on and 
wind-off conditions was the prior experience in the work conducted by NASA-Ames, the 
University of Washington, and McDonnell-Douglas. In each case, one hundred images 
was adopted in all work done to date. 

The ten averaged images of the wind-off condition were loaded into ten 
frame buffers of the EPIX framegrabber board. As with the initial averaging of the raw 
image data, the ten buffers were averaged together to form a single image. This single 
image, a result of averaging one hundred "raw" images together, was then saved to floppy 
disk. The corresponding wind-on images were averaged in a similar manner, and the 
resultant image saved with its wind-off counterpart. 

The procedure used in averaging the data in the present work was slightly 
different than used at NASA-Ames, the University of Washington, and McDonnell- 
Douglas. In those cases, one hundred images were collected and then averaged in one 
step. In the present work, ten averaged images were averaged to form the final image. 
Again, this is a result of the limited memory available on the frame grabber board. 
However, since each image has a random distribution of data associated with it, the 
averaging technique used will have little effect on the outcome of the result. 

The significance of image averaging can be seen in Figures 28 through 30. 


In each figure, an image and a graphical distribution along a single line of pixel values in 
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the image is shown. Figure 28 shows a single raw image of the wind-on condition at 80 
psig. Figure 29 shows the effect of averaging ten wind-on images together and Figure 30 
shows one hundred images averaged together. In each figure, an image and a graphical 
distribution along a single line of pixel values in the image are shown. 

(2) Dark Current. Appendix D details what dark current ıs and how it affects 
the storage of photogenerated carriers in the potential well of a photo-sensitive array 
element. The effect of the dark current on an image is that it produces an offset in which 
the zero incident light condition does not correspond to the zero grey level. To correct 
for this offset, a "dark" image was acquired and then subtracted from the wind-off and 
wind-on conditions. 

To acquire a dark image, a lens cover was placed on the lens prohibiting 
any light from striking the CCD array, and then an image was acquired. One hundred 
images were taken and averaged as was done with the wind-on and wind-off conditions. 
A dark image was generated after each of the wind-on runs. Analysis of the dark current 
image showed that 1 to 2 grey levels were present, which were the levels subsequently 
subtracted from the wind-off and wind-on images. 

(3) Registration. When a pixel-to-pixel comparison of two images of the 
same object field taken from the same sensor at different times is desired, it is necessary 
to spatially register the images to correct for relative translational shifts, rotational shifts, 


and geometrical distortions. In the general case, this is accomplished in the software 
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Figure 28a. Single Image of the Underexpanded Jet Across the Flat Plate at 80 psig 
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Figure 28b. Graphical Representation of a Single Line of Pixels of a Single Image 
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Figure 29a. Ten Averaged Images of the Underexpanded Jet Across the Flat Plate 
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Figure 29b. Graphical Representation of a Single Line of Pixels of Ten Averaged 
Images 
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Figure 30a. One Hundred Averaged Images of the Underexpanded Jet /Flat Plate 
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using a pair of transformation equations relating the new coordinates to the old 
coordinates. An example set of transformation equations, termed the affine projection, 
are 

X — (ay X + apy Гах та, у“) 
and 

у = (aX + apy +2,)Ma, x + ayy +1) 
where x y are the old coordinates and x',y' are the new coordinates. The eight unknown 
coefficients in the two equations can be determined if the coordinates of four points, called 
control points, can be determined from each of the two images involved. [Ref. 21:p. 71] 

In the jet/plate interaction experiment, the stiffness of the plate was 
expected to be sufficient to withstand the forces imposed on it by the jet so that any 
deformation that did occur would be insignificant, making registration unnecessary. 
However, comparison of the wind-on and wind-off images showed that the plate did 
experience a horizontal translation when the jet was on. It was found that the thrust 
produced by the jet caused a deformation in the structure supporting for jet apparatus. 
Thus, the jet apparatus and flat plate in combination experienced a displacement relative to 
the camera. 

Examination of the wind-off and the corresponding wind-on images 
showed that for the 65 and 80 psig cases , the images had shifted one pixel width in the 
direction of the thrust, and two pixel widths for the 90 psig case. The image-registration 


algorithm in the EPTX software was found to be unsuitable for this correction. In 
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addition, the lack of "natural" features on the plate to serve as control points did not allow 
the use of the software. [When the object lacks such features and image registration is 
anticipated, artificial marks such as black circles can be distributed on to the surface. 
In the present case, registration was not expected to be needed, the plate was not so 
equipped. | 

The registration between the wind-off and wind-on images was a 
relatively simple issue involving rigid body motion in one dimension. For this particular 
situation the images were aligned using the "Copy and Shift Image" function in the EPIX 
software which allowed an image to be shifted in its buffer up or down, left or right by an 
operator- specified number of pixels. The registration scheme was accomplished by first 
loading the wind-on and wind-off images into consecutive image buffers. The wind-on 
condition was considered the "warped" image and the wind-off image the "reference" 
image to which the warped image would be registered. A cursor overlay was 
simultaneously displayed on the image monitor with the wind-off image. The cursor was 
adjusted so that its horizontal character overlaid the pressure ports along the centerline of 
the plate while the vertical character of the cursor overlaid the pressure ports along the 
vertical centerline of the plate. The cursor remained fixed relative to the screen of the 
monitor; this served as a reference for comparison of object feature location when shifting 
between the wind-on and wind-off buffers. The wind-on image was then shifted relative 
to the vertical-cursor character until the vertical centerline pressure ports were overlaid by 


the vertical cursor character. When complete, the images were considered registered. 
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The accuracy of using this approach was judged to be similar to that of 
registration software since the registration software normally requires the operator to 
manually locate the control points in each image. In addition, since the body moved 
relative to the light source, consideration must be given to whether changes in brightness 
were due to the pressure variations, or were caused by the motion of the plate. 
McClachlan ef al [Ref. 22:p. 3] indicate that this is of greatest concern when model 
deflections and distortions are extreme. Since the flat plate motion consisted of rigid body 
translation, the magnitude of which was extremely small, this concern was considered to 
be insignificant. 

(4) Ratioing. While the reason for ratioing the wind-off image with the 
wind-on image was to satisfy the Stern-Volmer relation and produce a pressure map, there 
were other benefits which arose from the procedure. In Figure 31 it is easy to discern the 
gross features of the pressure distribution in the wind-on image. The ratioing process had 
the affect of enhancing the finer detail in the image that would otherwise be too subtle for 
detection. [Ref. 7:p. 7] 

Examination of the wind-off image in Figure 32 shows the evident 
brightness variations from point to point across the surface. This was due to a 
combination of differences in illumination, paint thickness variations, and features in the 
model itself. Figure 33 shows a three dimensional graphical representation of the wind-off 
intensity values across the surface of the plate. This illustrates the variations that exist 


across the plate. The result of ratioing was the elimination of these non-uniformities. 
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Figure 31. Wind-on Image of the Underexpanded Jet at 80 psig Across a Flat Plate 
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Figure 32. Wind-off Image for the 80 psig Condition 


65 





Т е ШЫ, 
i ETTR мж 


MR Рр 


CEDAR 
oL 





Figure 33. Three Dimensional Intensity Plot of the Wind-off Condition of the Plate 
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А The "Ratio Image Operation" under the "Two Image Arithmetic" menu of 
the EPIX software formed the weighted ratio of corresponding pixels. The operation is 
mathematically expressed as 

(с0 * PıxB + c1)/(c2 * PıxA + c3) 
where the coefficients cO, c1, c2, and c3 are modifiable. For the current work, PixA was 
designated the wind-on image while PixB, the wind-off image. In addition, the coefficient 
cO was assigned the value of 100 while the coefficients c1, c2, and c3 were assigned the 
values 0, 1, O respectively. The reason for assigning cO the value of 100 was 
to scale the resultant ratioed intensity map such that differences of the grey values between 
neighboring pixels could be discerned. 

Figure 34 illustrates the ratioed image for the underexpanded jet across 
the flat plate at 65 psig plenum pressure. The importance of proper registration of the 
wind-on and wind-of images prior to ratioing cannot be overstated. To illustrate the effect 
of improper registration, the images used to form the resultant image in Figure 34 were 
modified where the wind-on condition was shifted an additional five pixel widths in the 
direction of the thrust. The result of ratioing these two images is shown in Figure 35. A 


reduction in fine details, as compared to Figure 34, is clearly evident. 
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Figure 34. Ratioed Image of the Wind-off and Wind-on Images for the Flat Plate 
/Jet Interaction at 80 psig 
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Figure 35. Result of a Ratioed Image Following Improper Registration Technique 
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b. Calibration 

Two methods were available to calibrate the intensity ratio map, the intent of 
which was to obtain the sensitivity coefficients for the Stern-Volmer equation and 
establish the pressure map. 

The first method, termed "a prior", involved calibrating a paint sample in a 
controlled pressure and temperature chamber. Using the data obtained from the chamber, 
and knowing the temperature of the model surface, the appropriate values for the 
Stern-Volmer coefficients could be found [Ref. 7:p. 7]. 

The second method, termed "ir situ", which was the method employed in the 
present work, used the pressure obtained using a conventional pressure tap in the model 
surface and compared this with the intensity from the luminescent paint adjacent to the 
tap. Using data from a number of taps, the sensitivity coefficients were obtained by a 
least-squares linear fit. The coefficients A and B, being functions only of temperature, 
were good over that portion of the surface which was at the same temperature as the area 
containing the pressure taps. [Ref. 7:p. 10] 

The approach used to produce pressure data from the ratioed image was 
similar to that used by Kavandi [Ref. 9:p. 94]. This entailed averaging the data from five 
adjacent rows of pixels near the pressure taps to produce a single row of data. (This, has 
the effect of averaging a total of five hundred images together.) 

To do this, the ratioed image was loaded five times into consecutive frame 


buffers of the frame-grabber board. Each successive image was then vertically offset by 
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one pixel row with respect to the previous image, 1.e., the first image remained 
unmodified, the second image was shifted one pixel row down in its frame buffer with 
respect to the first image, the third image shifted two pixel rows down with respect to the 
first image, etc. The five buffers were then averaged. Each row of pixels in the resultant 
image, with the exception of the first five and the last five rows of the image, was then an 
average of the corresponding row in each of the five images. A row of pixels close to the 
pressure taps was then selected to serve as calibration data. Figure 36 shows the data 
from the selected line of pixels for the flat plate in the underexpanded jet at 80 psig, prior 
to averaging with four adjacent lines of pixels. Figure 37 shows the data for the same line 
after the averaging procedure. [It is noted that the "spike" on the left side of the graph is 
due to paint that has separated from the leading edge of the plate]. From an examination 
of Figure 36 and Figure 37 it can be concluded that the five-row averaging procedure had 
the effect of reducing the noise level in the distribution with negligible effect on the 
distribution itself. 

The calibration curves that resulted from plotting intensity values against the 
measured surface pressure at jet stagnation pressures of 65, 80 and 90 psig, are shown in 
Figures 38, 39 and 40, respectively. All data were obtained using the line of pixels and the 
line of pressure taps along the centerline of the jet. Two separate linear fits are shown in 
each of the figures. This was done when it was realized that the data were poorly 
represented by a single curve, and that large variations in temperature over the surface of 


the plate was the probable cause. It is possible that the groups of points which correlate 
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Figure 36. Calibration Line of Pixels for the 80 psig Case Prior to Spatial Averaging 
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Figure 37. Calibration Line of Pixels for the 80 psig Case After Spatial Averaging 
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Callbration Curve For The Underexpanded 
Jet Nozzle Across A Flat Plate 
(65 psig Plenum Pressure) 
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Figure 38. Calibration Curve for the Underexpanded Jet/Flat Plate at 65 psig 
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Calibration Curve For The Underexpanded 
Sonic Jet Across A Flat Plate _ 
(80 psig Pienum Pressure) 
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Figure 39. Calibration Curve for the Underexpanded Jet/Flat Plate at 80 psig 
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Calibration Curve For The Underexpanded 
Sonic Jet Across A Flat Plate 
(90 psig Pienum Pressure) 
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Figure 40. Calibration Curve for the Underexpanded Jet/Flat Plate at 90 psig 
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linearly are from regions of similar surface temperatures. However, the complexity of the 
flow over the plate surface, with a three-dimensional separation contained within it, makes 
the estimation of surface-temperature distribution very difficult. These results clearly 
indicate the need to measure surface temperature if in-situ calibration 1s to be used. 

c. Field Pressure Map 

Pseudo-coloring (or "false coloring") of a gray scale image was used to 

accentuate areas of an image that are not easily discerned by the eye. This was 
accomplished by providing an assignment of amounts of red, green and blue to a 
corresponding value of intensity in the gray-scale image. Through the use of look-up 
tables, the software allowed specification of linear ramps, logarithmic curves, polynomials, 
gamma correction, inversion and transposition in a pattern or sequence of a non-specific 
nature. For a gray-scale image the index value is synonymous with its intensity value. 
Reference 23 may be consulted for further information concerning the manipulation of 
lookup tables. The color map or sequence that was developed for the present work is 
shown in Figure 41. The constants A and B represent the minimum and maximum index 
values of the interval over which the color map was defined. All other index values were 
assigned the pixel value of zero in the color map. The Бий and maximum index 
values were found by conducting a histogram on the gray-scale image which gave a 
breakdown of the pixel intensity values and the frequency of occurrence in the image. The 
intensity values associated with the pressure-sensitive paint images are usually grouped in 


an interval between the 0 to 255 index value range. 
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Figure 41. Graphical Representation of the Red, Green, and Blue Lookup Tables 
The constant C represents the mid-point of this interval. The constants D and E were 
varied to obtain the desired bandwidth of the green lookup-table input which influenced 
the mid portion of the color spectrum. 
The final pressure maps for the 65, 80, and 90 psig plenum conditions are 
given in Figures 42, 43, and 44 respectively. Although quantitative information is not 
available, the results in Figures 42, 43, and 44 show good qualitative agreement with their 


schlieren counterparts in Figures 12, 13, and 15 respectively. 
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Figure 42. Pressure Map of the Underexpanded Jet at 65 psig Across a Flat Plate 
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Pressure Map of the Underexpanded Jet at 80 ps 
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Pressure Map of the Underexpanded Jet at 90 ps 
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IV. APPLICATION TO SHOCK-BOUNDARY LAYER INTERACTION IN A 


SUPERSONIC WIND TUNNEL 


A. DESCRIPTION OF FACILITY 

The supersonic wind tunnel, located in Bldg. 216 of the Gas Dynamics Laboratory 
[Bldg. 216] at the Naval Postgraduate School, was of the blow-down type. Supply air 
was produced by the same supporting facility as the underexpanded sonic jet apparatus. Á 
schematic of the tunnel and supporting facility is shown in Figure 45 while Figure 46 
shows a photograph of the wind tunnel. 

Air to the tunnel, from a maximum pressure of 300 psi, was controlled by a pneumatic 
control valve. The control valve used a diaphragm actuator. Regulated instrument air 
from the supply system was routed to one side of the diaphragm and the controlled 
pressure, from a tap on the plenum, was connected to the other. This allowed precise 
control of the plenum pressure. A pressure gauge installed on the plenum provided 
pressure readout. 

The plenum chamber was cylindrical with an internal diameter of approximately 20 
inches. With the exception of a splash plate mounted perpendicular to the flow near the 
entrance of the chamber, the plenum was empty. The flow exited the chamber through a 


contraction and circular-to-rectangular transition section, into the convergent-divergent 
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Figure 46. Supersonic Wind Tunnel Apparatus 
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nozzle. The nozzle and test section,shown schematically in Figure 47, were comprised of 
a set of two (interchangeable) aluminum blocks forming the contour and two flat plates 
forming the side walls and enclosing the constant-width section. 

The test section was nominally four inches in width and four inches in height. Six-inch 
diameter circular glass windows on each side of the test section provided optical access. 
Two small gate valves were mounted on each of the side walls just downstream of the test 
section. They were used to allow "bleed" and thus control back pressure in order to 
position the "starting" normal shock in the test section. The flow exited the test section 


into an exhaust duct which vented to the atmosphere outside the building. 
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Figure 47. Nozzle Blocks and Test Section Exploded View 
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1. Nozzle Blocks 


The interchangeable pairs of nozzle blocks were each designed to produce a specific 
Mach number in the test section. For the present study, two Mach numbers, 1.4 and 1.7, 
wer used, since this was the range of Mach number of interest in the proposed study of 
control of shock-boundary layer interaction. A set of blocks for M=1.4 existed. 
However, 1t was necessary to design and manufacture a set of blocks to provide M=1.7. 
The program described in Appendix E was used to derive the contour for the M=1.7 
nozzle. The new blocks were otherwise similar to the existing blocks and were 
manufactured to the orıginal dramings using the new contour. 
a. Performance Evaluation of the Mach 1.7 Nozze Blocks 

To evaluate the new nozzle blocks, several runs at varying plenum pressures 
were conducted. A series of pressure taps along one side of the tunnel allowed pressure 
data to be collected. The taps extended approximately four inches upstream of the nozzle 
throat and downstream to the entrance of the test section. A metal blank containing 
additional taps was installed in place of one of the test section windows to measure 
pressures through the test section. One tap downstream of the test section measured the 
pressure of the exiting flow. 

The data acquired indicated that the normal shock associated with the nozzle 
starting process rapidly passed through the test section into the exhaust duct as soon as ' 
the operating pressure ratio decreased below the first critical condition. Unlike with the 


M=1.4 blocks, the bleed valves on the sides of the tunnel could not be used to position the 
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shock. It was realized that some form of blockage would have to be introduced into the 
passage at some point downstream of the test section in order to "set" the shock in the test 
section. 

The mounting points for the gate valves, which were located three inches 
downstream of the test section and were aligned such that their centerlines coincided with 
each other, were adapted so that circular rods of varying diameter could be installed 
across the tunnel. A total of four rods was used. Operating the tunnel at 75 psi plenum 
pressure, the static pressure distribution was recorded for each configuration. The effect 
of rod diameter on shock position is clearly seen in Figure 48. 

The 3/8 inch diameter rod was found to place the shock just prior to the test 
section. To position and control the shock precisely in the center of the test section, a 
means of bleeding off mass flow was developed. The solid rod was replaced by a 3/8 inch 
outside diameter tube. A series of 11/64 inch holes was drilled through the tube wall, in a 
line along its 4-inch length. The tube was then mounted between the gate valves with the 
holes aligned with the flow. This arrangement allowed a controlled mass flow to pass 
through the holes into the tube and out of the tunnel via the gate valves. The gate valves 
could be adjusted to vary the amount of outgoing mass, thus moving the shock to the 
desired position in the test section. A schematic of this arrangement is shown in Figure 49 


while Figure 50 shows a photograph of the bleed-off tube mounted in the tunnel. 
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Figure 48. Tunnel Pressure Distributions as a Result of Varying Rod Diameters 
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Figure 49. Cross-Sectional View of the Mass Bleed-Off System 
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Figure 50. Mass Bleed-Off Tube Mounted in the Tunnel 
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B. INSTRUMENTATION 

The metal blank, containing eleven pressure taps, and used as described above to 
measure the pressure distribution in the test section, was used as the model surface on 
which to observe the interaction between the starting normal shock and the nozzle 
boundary layer using pressure-sensitive paint. The metal blank was installed on one side 
of the test section and a six inch glass window was installed on the other. Concern for the 
spectral transmittance of the glass window led to an analysis of the glass for the range of 
350 nm to 700 nm using a spectrometer. The results of the test, shown in Figure 51, 
indicated that the transmittance levels were satisfactory for the current work. 

The optical measurement system was the same as was used for the underexpanded jet 


study described in chapter three. A schematic of the setup is shown in Figure 52. 


C. EXPERIMENTAL PROCEDURE 

It was found that the shock-boundary layer interaction in the test section was 
inherently unsteady. This precluded the use of the procedure outlined in chapter three for 
acquiring the wind-on data. That method required the event of interest to be steady and 
repeatable, thus enabling one hundred images to be taken over a series of ten runs. The 
approach taken for the unsteady process was simply to create as many image buffers as 


possible for a single sequence, and to execute this sequence in the shortest period of time. 
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Figure 51. Spectrometer Results for the Optical Window 
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Pressure Measurement set-up Shock-Boundary Layer Interaction Studies 
Gas Dynamics Laboratory 
Figure 52. Instrumentation Setup for the Supersonic Wind Tunnel, Shock- 
Boundary Layer Experiment 





The maximum number of image buffers that could be created for image capture was a 
function of video resolution. To maximize the number of available buffers it was 
necessary to minimize the video resolution. To accomplish this, the camera was focused 
on the test section so that the surface of the metal blank was clearly displayed on the video 
monitor. The video resolution was set at maximum (752 x 480) at this point. With the 
video resolution menu displayed on the computer screen, the parameters controlling the 
numbers of pixels and pixel lines were iteratively decreased until the camera's field of view 
shown on the video monitor exhibited only the test surface, the rest of the field of view 
being deleted. The resulting resolution was 600 pixels per line and 340 lines of video or 
600 x 340 (assuming "Interlace: Use Distinct (2) Fields?" is active) which allowed 21 


image buffers to be available for image capture. Since McDonnell-Douglas [Ref. 6:p. 8] 
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had shown that averaging data from as few as eight images significantly reduced the noise 
on the signal, it was concluded the 21 images would be sufficient. 

The period required for a sequence to be completed was based on the time interval 
specified between consecutive samples, the shortest interval being 1/30 of a second. Since 
a video frame was generated every 1/30th of a second by the camera, the total time to 
execute a sequence containing 21 image buffers was 1.37 seconds. 

The experiment began with the Mach 1.4 blocks installed in the tunnel. Using the 
shadowgraph technique, the normal shock was positioned in the test section by adjusting 
the backpressure through the gate valves, while the plenum pressure was set at 30 psi. 
When the valves were properly adjusted the tunnel was shut down, the metal window 
blank with pressure sensitive paint applied was installed, and the camera and illumination 
system were located eight inches away from the optical glass window. A real-time image 
of the test section was displayed on the image monitor, and the parameters to execute the 
desired sequence described above, were input into the computer. A cover was placed 
over the test section to prevent photodegradation of the paint, the light source was turned 
on and allowed to "ramp up" and stabilize. The wind tunnel was then started, the plenum 
pressure set at 30 psi, and the laboratory lights were turned off. The wind-on data were 
captured while pressure data were taken simultaneously using the CALS YS2000 data 
acquisition system. Twenty samples per port were taken at a sampling rate of 10,000 Hz. 
Three such runs were conducted. Between each run, wind-off and dark-current images 


were collected. Temperatures for the tunnel and the test plate were not available. 
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The Mach 1.7 nozzle blocks were then installed with the mass bleed-off system 
described earlier. With the plenum pressure set at 70 psi the shadowgraph was again used 
to position the normal shock. With the test plate installed, three runs were conducted in a 
similar manner to that described for the M=1.4 blocks. The data were post-processed 


using the techniques described in chapter three. 


D. DISCUSSION AND RESULTS 

The flow structure expected to be produced by the interaction between a normal shock 
and a turbulent boundary layer is shown schematically in Figure 53. The 
pressure-sensitive paint technique was to be used to measure the surface pressure through 
the interaction region. Previous work by Perretta [Ref. 24] involved a study of the flow 
across the shock structure shown in Figure 53, and the resulting separation, for the given 
M-1.4 tunnel configuration used in the present study. 

For the M=1.4 case, the pressure map obtained is shown in Figure 54. It can be seen 
that the interaction was nearly two-dimensional. To illustrate the distribution of the 
pressure rise, a line of pixel intensity values taken along the plate is shown in Figure 55. 
The profile resembles those of Kooi [Ref. 25:p. 30-8] for turbulent boundary layer 
separation on a flat plate due to normal shock interaction. The steepest slope of the curve 
indicates the start of the interaction, according to Kooi's work. The separation point can 


be identified as that point in the curvature where either a discontinuity in the curvature 
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Figure 53. Schematic of the normal shock-boundary layer interaction 
exists or the degree of change in curvature is much less compared to that of the start of 
the interaction, indicating a much slower rise in pressure. From the distribution of the 
intensity values, this point 1s difficult to determine. This can be attributed to oscillations of 
the normal shock, and the pressure-sensitive paint's poor dynamic response to these 
changes, causing an averaging effect and eliminating changes in curvature from the curve. 
Also, the time to execute the data collection sequence was much longer than the 
oscillation period of the shock, adding to the error in the results. The absence of a plateau 
in the intensity distribution suggests that the separation region was very short. 

In the case of M=1.7, the determination of separation is more definite. The pressure 


map is shown in Figure 56 while the corresponding streamwise line of pixel-intensity 
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Figure 54. Mach 1.4 Pressure Map of the Shock-Boundary Layer Interaction 
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Figure 55. Graphical Representation of the Pressure Distribution of the Shock- 
Boundary Layer Interaction at Mach 1.4 





values is displayed in Figure 57. The intensity curve here shows a definite plateau, a 
region of nearly constant pressure, which is an indication of an extended length of 
separation. This region corresponds to the area illustrated by the narrow green band in 
the pressure map of Figure 56. 

The calibration curves for the M=1.4 and M=1.7 data are given in Figures 58 and 59 
respectively. In contrast to the calibration curves developed for the -—" jet/ 
flat plate data, areas of similar temperature are more easily identified in the 
shock-boundary layer experiment. The data points upstream and downstream of the 
normal shock have been linearly approximated separately since there was a near 
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Figure 56. Mach 1.7 Pressure Map of the Shock-Boundary Layer Interaction 


98 





Value 


эх 
w 
х 
рю 
де 


% Coordinate 
Figure 57. Graphical Representation ofthe Pressure Distribution of the Shock- 
Boundary Layer Interaction at Mach 1.7 





discontinuity in temperature and heat transfer at the shock impingement location. There is 
seen to be a change in the sensitivity of the paint to pressure change which occurs at the 
shock. The paint is more sensitive downstream of the shock, which is consistent with the 


occurrence of higher temperature. 
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Calibration Curve For The Shock-Boundary 
Layer Interaction In The Supersonic Windtunnel 
(Mach 1.4 -- 30 psig Plenum Pressure) 
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Figure 58. Calibration Curve for the Mach 1.4 Interaction 
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Calibration Curve For The Shock-Boundary 
Layer Interaction In The Supersonic Windtunnel | 
(Mach 1.7 -- 75 psig Plenum Pressure) 
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Figure 59. Calibration Curve for the Mach 1.7 Interaction 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The pressure-sensitive paint (PSP) technique offers a flexible and convenient method 
for measuring surface static pressures in aerodynamic testing. Its most significant feature 
is the non-intrusive manner in which it provides a mapping of a continuous pressure field 
over a surface. Because of its detailed spatial resolution (when compared to conventional 
pressure taps), a complete picture of the events occurring on the surface is available, 
including details not necessarily anticipated prior to testing. It appears to be highly 
suitable as a technique to use to validate results calculated using computational fluid 
dynamics (CFD). In the present study, pressure-sensitive paint was applied successfully in 
two experimental configurations involving shock-boundary layer interaction. The 
conclusions drawn from the present application are relevant to any first-time application of 
the technique. 

The luminescent technique, although conceptionally simple, required learning. It was 
important to understand the entire process, including each of the components comprising 
the pressure measurement system (paint, optics, camera, software), the limitations of the 
system (particularly of the luminescent paint), and the interpretation of the measurements 
through calibration, before meaningful data were acquired. Several attempts to progress 
through the entire process were found to be necessary. Thus, the initial experimental 
arrangement (a flat plate immersed in an underexpanded jet) served well as a test bed for 


building experience before more complicated experimental geometries were attempted. 
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The components that comprised the current measurement system were adequate for 
initial experience, system development, and measurements. However, in order to make 
highly quantitative measurements routinely, some refinements are recommended. The 
camera and the digitizer combination are, by far, the most critical components in the PSP 
measurement system. The accuracy and precision of the measurements is highly 
dependent upon the dynamic range and signal-to-noise ratio of these components. The 
camera and digitizer used had a grey-scale resolution of eight bits giving a dynamic range 
of 256:1. Since many measurements, such as in shock-boundary layer interaction 
experiments, may generate changes in the luminescent intensity as little as two percent, it 
is crucial that the camera and digitizer be sensitive enough to these changes and to assign a 
gray-scale value that 1s appropriate for the given intensity level. Recommendations from 
NASA Ames Research Center and McDonnell Douglas have suggested a minimum of 14 
bits for gray level resolution. 

With a relatively low signal being detected, the signal-to-noise ratio at high light levels 
is limited by the photon shot noise. The accuracy of measurement is contingent upon the 
highest possible full-well potential of the CCD element, which requires a minimum 
generation of dark current. Active cooling of the CCD array would maximize the 
signal-to-noise ratio of the array. Such developments in CCD cameras can be expected 
[For example, Princeton Instruments Inc. currently offers a Penta Max camera that uses 


air circulation to cool its CCD array, which achieves -45 degrees Celsius]. Similarly, 
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system noise can be reduced by combining the digitizer with the camera, thereby reducing 
electrical noise, and then feeding the signal directly into a memory buffer in the computer. 

The EPIX 4Meg Model 12 frame grabber and memory buffer and its associated 
software was limited in data collection. Since the board was limited to capturing only 11 
images per sequence (assuming maximum resolution), the acquisition of one hundred 
images was far to slow to avoid inaccuracy of the results due to photodegradation of the 
paint. The EPIX 48Meg Model 12 frame grabber, for example, would provide the ability 
to capture up to 132 images per sequence. Not only would this reduce or eliminate the 
error associated with photodegradation, but image processing time would be significantly 
reduced. Data acquisition speed is mandatory in experiments in which the flow 15 not 
completely steady. 

A similar issue concerns the present image-processing software. Designed for general 
applications, it is not the most efficient, or the most suitable, for the PSP applications. 
Although single-function algorithms such as registration software to spatially align two 
images together, to compensate for temperature sensitivity and apply the calibration to the 
data, and routines to pseudo-color the results must be written, consideration should be 
given to automating the entire process, from acquisition of the data to production of the 
calibrated pressure field map. 

The temperature sensitivity of the paint and its effect on the results is a major concern 
if pressure measurements are to be made in the compressible flow regime. To accurately 


measure the pressures on a surface, it will be necessary to know the surface distribution of 
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temperature so that the calibration curves can be adjusted accordingly. Although 
temperature measurements could be made using thermocouples, this approach raises the 
same issues as those associated with using pressure taps to make pressure measurements. 
A possible solution is to use, concurrently, a high-performance infrared camera to record a 
high-resolution, continuous image of the temperature distribution on the surface. Since 
the results from infra-red cameras are produced using a digital process similar to that used 
in PSP, it is possible that the techniques could be combined to produce a temperature- 
compensated pressure map. 

Finally, the applicability ofthe present technique to more complicated model 
configurations with more limited optical accessibility of the surface of interest, has yet to 
be attempted here. The combination of curved surfaces with significant temperature 
variations, and with limited optical access, will provide a significant challenge to the 


successful quantitative use of PSP. 
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APPENDIX A. RUSSELL-SANDERS COUPLING 


Associated with the angular momentum of an electron 1s a magnetic moment. When 
the electron is placed in a magnetic field, the interaction of the field with the magnetic 
moment of the electron causes the angular-momentum vector to precess about the 
direction of the external field, much the same way a toy top precesses in a gravitational 
field. However, in contrast to the top, which can have an infinite number of momentum- 
vector orientations with respect to the field direction, the angular-momentum vectors of 
the electrons in an atom (or molecule) are allowed only certain orientations. The solutions 
of the Shródinger equation dictate permitted orientations; not every orientation is allowed. 
The only allowed directions are those at which the components of the angular momentum 
along the direction defined by the magnetic field have certain quantized values. This 
behavior is called space quantization. [Ref. 10:p. 491] 

Each orientation of the momentum vector with respect to the external field direction 
corresponds to a unique energy state. Since the angular motion of the electron is 
equivalent to an electromagnet, the energy of the interaction between the magnetic 
moment it possess and an external field can be expressed as E = u H cos 0, where H is the 
strength of the external field, O is the angle between the field and momentum vector, and 
H is the magnetic moment. [Ref. 26:p. 25] 

Now consider what happens when a single unpaired electron is not in the presence of a 


magnetic field. Here, the orbital angular momentum of the electron creates the magnetic 
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field, defining a specific direction in space. Thus the angular momentum associated with 
the intrinsic spin of the electron is oriented with respect to this magnetic field. According 
to space quantization principles, the spin angular-momentum vector of the electron can 
have two orientations with respect to the direction defined by the magnetic field; this 
corresponds to the doublet line structure observed for the spectra of an atom, and is the 
basis for the spin quantum number (+1/2 and -1/2) assigned to the electron. [Ref. 26:p.27] 
In an atom with several electrons, the orbits are usually oriented first with respect to 
each other, giving a resultant orbital magnetic moment. Similarly, there is a resultant spin 
moment which is oriented with respect to the resultant magnetic field due to the orbital 
motion. With space quantization, the resultant spin vector may have several orientations, 
each corresponding to a different energy level. Thus, the atom (or molecule) may have 
several energy levels for a single configuration. [Ref. 27:p. 79] These multiple energy 
levels for a single configuration of an atom (or molecule) is the basis for Russell-Sanders 


coupling which impacts the multiplicity of state that an atom or molecule can possess. 
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APPENDIX B. LUMINESCENCE COATING CHARACTERISTICS 


The luminescent coating used in the present work possessed several undesirable 
characteristics of which the user must be aware if accurate results are to be obtained. A 
summary of these is presented below. For a more in-depth discussion, the reader may 


refer to the references mentioned. 


A. TEMPERATURE SENSITIVITY 


The photoluminescence process is temperature sensitive by nature. The coating used 
in the present work exhibited a non-linear response, and intensity decreased with 
increasing temperature. This response at a constant pressure is shown graphically in 
Figure B1. [Ref. 1:p36] The temperature dependence has a significant effect on the 
calibration curves of intensity vs. pressure. As is shown in Figure B2a, the variations in 
slope can be significant [Ref. 1:p. 62]. 

An interesting property displayed by the calibration curve for the coating is the 
condition where the coating-sensitivity coefficients add to unity (i.e. A + B = 1, in the 
Stern-Volmer equation). This occurs when the data taken for I and I, are taken at the 
same temperature. If this is not the case, the coefficients add to other than one, 
A+B#1. The consequences are illustrated in Figure B2. In Figure B2(a) only the 23.7 ° 
C curve exhibits the property that A + B = 1 while the 50 and 6 °C calibration curves were 


calculated with I and I, at different temperatures. In Figure B2(b), the 50 and 6° C 
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Figure B1. Luminescence Intensity as a Function of Temperature [From Ref. 1:p. 36] 
calibration curves have been replotted using data for both I and I, obtained at similar 
temperatures. Clearly, the slope and zero intercept have been significantly altered. [Ref. 


1:p.37] 


B. PHOTODEGRADATION 
When subjected to ultraviolet light, the luminescent coating undergoes 


photodegradation, resulting in a decrease of its response with time of exposure. A 
representative example of this is shown in Figure B3 where intensity is plotted as a 
function of time of exposure. In this case the coating lost 46% of its response after one 


hour of exposure. [Ref. 1:p. 36 ] 
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Figure B2. Effect of Temperature on Coating Calibration Curves [From Ref. 1p: 37] 
(a). I, for Curve Measured at 23.7” C 
(b). I, for Curve Taken at its Respective Temperature 
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Figure B3. Illustrative Example of the Effect of Photodegradation on 
Response Time [From Ref. 1p:36] 
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C. TIME RESPONSE 

The luminescent coating used in the present work has been found to exhibit a time 
response on the order of 1 second limiting measurements of pressure fluctuations to no 
greater than 1 Hz. [Ref. 8:p. 3] 

This limitation is not a characteristic of the active molecule. In theory, the time 
response of the coating to pressure transients may be as fast as the decay of luminescence. 
For platinum octaethylporphyrin this is approximately 0.1 milliseconds. The problem lies 
with the polymer matrix in which the active molecule is embedded. By limiting the rate of 


oxygen diffusion, the polymer slows the optical response to surface pressure fluctuations. 


[Ref. 1:p. 41] 


D. COATING THICKNESS / INDUCTION PERIOD 

Kavandi [Ref. 9:p.86] has conducted tests to determine 1f film thickness has any effect 
on calibration results. Over the range of 5.5 micrometers to 17 micrometers, consistent 
with normal application, no can variations were noted. However, a noticeable lag, 
or "induction period", existed between the time excitation light was turned on and the time 
that the detected luminescence reached maximum intensity. The period of maximum 
intensity was observed to be in the range of 30 to 45 seconds and it is therefore advisable 
to allow this period to elasped before taking data. A graph of absolute intensity as a 
function of time following exposure to the exciting light is shown in Figure BA. [Ref. 


9:p.89] 
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Figure B4. Induction Period of the Luminescent Coating with Thicknesses in 
Excess of 17 Micrometers [From Ref. 9:p. 89] 
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APPENDIX C. EPIX SOFTWARE SUMMA RY 


This appendix outlines the software steps used in acquiring, processing and reducing 
image data. All information discussed here may be found in [Ref. 23]. The appendix 1s 
intended to serve as a user guide but is not meant to be a substitute for the EPIX user's 
manual. The user's manual is vital to developing a sound understanding of the software in 
order to be able to use it to its fullest extent. The data acquisition and image processing 
system is set up such that the menus and commands are displayed on the computer 


monitor while the images are displayed on the Sony image monitor. 


A. IMAGE DATA ACQUISITION 

The main menu for the EPIX 4MIP-4MIPTOOL software is shown in Figure C1. 
Prior to data collection, the video resolution is set. Selecting "Video Resolution" from the 
main menu, the video resolution sub-menu will be displayed as in Figure C2. Check to 
insure that the "Interlace: Use Distinct (2) Fields ?" is set to "ON". The sampling 
resolution is considered next and it is application dependent. The sampling resolution sets 
the size of the video window displayed on the image monitor, allowing the selection of a 
smaller field of view. The maximum resolution 1s 752 pixels per line by 480 lines per 
frame or image (Two fields make up one frame; however the number of lines selected is 
based on each field, the maximum number being 240. For both fields to be active, the 


interlace function must be on). Upon setting the desired sampling resolution, the 


115 


4MIP 4MEG VIDEO Image Acquisition, Display, Processing, Analysis 
v2.8 Copyright (C) 1984-1993 EPIX, Inc. All rights reserved 


I! COMMAND >> |< Quit Menu 


Video Formats ? Help Key Usage 
Video Resolution > Image Test Patterns & Sequences 


Video Digitize/Display Image Processing 
Motion Sequence Capture/Display Image Printing 
Special Operations & Modes Image File Load/Save 


Contrast & Lookup Tables Image Measurements 
Pixel Peek, Poke & Plot Paint, Draw & Text Overlay 
Image Zoom & Pan Feature Finders 


MIPX Scripts DOS Escape 
Custom Menu Obscure Menus 





Figure C1. "Main Menu" for EPIX Interactive Image Analysis Software 


> Video Resolution 
I COMMAND >> |< Quit Menu 
Interlace: Use Distinct (2) Fields ? ON 


|! Horz: Pixels Sampled per Line 752 |! Vert: Lines Sampled per Field 240 
|! Horz: Pixel Width 1 |! Vert: Line Height 
i! Horz: Left Edge 8 |! Vert: Top Edge 
і“ Horz: Right Edge 761 Vert: Bottom Edge 
i! Horz: Set L/R Centering ON |! Vert: Set Top/Bot Centering 
I Horz: Set Max Samples per Line OFF |! Vert: Set Max Samples per Field 


Variable Sampling (Pixel Width or Line Height) Parameters: 
|!Horz: Width Multiple at Edge 1 |!Vert: Height Multiple at Edge 
I!Horz: Width Unaffected Center 1 |!Vert: Height Unaffected Center 


¡Split Digitize/Display Screen Parameters: 
|! Horz: Division Position 2 |! Vert: Division Position 
! Horz: Division Spacing 8 |! Vert: Division Spacing 


Image Buffers in Memory 





Figure C2. "Video Resolution" Sub-Menu 
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number of frame buffers available for capturing images 1s displayed in the lower night 
corner of the menu. For 752 x 480 resolution, 11 frame buffers are available. After 
setting the resolution, return to the main menu by selecting "Quit Menu". 

From the main menu, select the "Video Digitize/Display"; this menu is shown in Figure 
C3. Select the function "Digitize", the result of which ıs a real-time video display of the 
camera's field of view. Adjustments to the camera can be performed, such as focus, gain, 
iris setting, etc., to achieve the desired image. When complete, return to the main menu; 
the real-time image will remain displayed on the monitor. 

Next select the "Motion Sequence Capture/Display" sub-menu; this menu is illustrated 
in Figure C4. The first step here is to choose the length of the sequence, and which 


buffers are to be used, by selecting the sequence starting and ending buffers. For example, 


I! COMMAND >> |< Quit Menu 


|i! Toggle Current Mode 

1! Display | ! Current Image Buffer 
! Digitize ! Image Memory Bit Write Mask 
! Blank Video ! Camera Control Mode 
! Stop Video Timing 


|! Set Genlock in Display ^ Switch Video after Odd Field 
i Set Genlock in Digitize 7 Switch Video after Even Field 
” Switch Video after Any Field 


Split Screen Modes 





Figure C3. "Video Digitize/Display" Menu 
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> Motion Sequence Capture/Display 


I COMMAND >> |< Quit Menu 


~ Sequence Starting Buffer 1 |” Available Image Buffers 
Sequence Ending Buffer 11 |! Current Image Buffer 
Trigger Options 


|! Digitize Sequence. Period: 
|! Display Sequence. Period: 
|! Digitize w. Sample Redisplay Sequence. Period: 
|! Digitize Continuous, Circular Sequence. Period: 


|! Reorder Circular Sequence 
! Average Image Sequence 





Figure C4. "Motion Sequence Capture/Display" Sub-Menu. 
in the flat plate/jet interaction work only 10 of the 11 buffers were used; thus the starting 
buffer was set to 1 and the ending buffer to 10. 

The time interval between consecutive samples is then set with the "Digitize Sequence. 
Period:" function. The period is based on units of video fields if the video is not interlaced 
(1/60 sec) or video frames if the video is interlaced (1/30 sec). The COHU camera was 
set internally to send an interlaced image to the frame-grabber board. Since the incoming 
signal to the board is interlaced, the software can be programmed using the interlace 
function to store either an image with both fields of the image (interlaced) or an image 
with a single field (non-interlaced). With the "Interlace: Use Distinct (2) Fields ?" set to 
active, the period will be based on units of frames. For the jet/plate interaction experiment 


collection, the period was set to the default value of 4 frame periods between each 
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captured image or 4/30 = 2/15 of a second between each captured image. For the 
supersonic wind tunnel experiment the period was set to 1 or 1/30 of a second between 
each captured image. 

One note concerning this function; to enter the desired period, the mouse or arrow 
keys are first used to position the cursor on the function. The period is entered using the 
numeric keys, the program responding to the entry by "highlighting" it red. The entry is 
then entered into memory by either strobing the left mouse button or striking the "enter" 
key. This action, however, also starts the execution of the algorithm for sequence image 
capture. When complete, all designated buffers will contain a captured image. If the 
image buffer used to display the real-time image lies between the starting and ending 
sequence buffers, it too will contain a captured image. Therefore, until the desired event 
is being displayed on the monitor in the real-time format, the entry for the period should 
not be strobed or the above sequence of events will have to be repeated to redisplay the 
real time image. Itis, of course, not necessary to display a real-time image to capture a 
sequence of images. By simply re-strobing the "Digitize Sequence. Period:" a new set of 
images is stored in the designated buffers. Displaying a real-time image of the desired 


event on the monitor is simply a matter of choice. 
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B. IMAGE PROCESSING 
1. Averaging 


The procedure outlined in Chapter III for data collection of the wind-on, wind-off, 
and dark-current conditions involved executing the above-described steps ten times, 1.e., 
ten images per sequence, executing the motion-sequence algorithm ten times, for a total of 
one hundred images. After each sequence run, the images in the designated buffers were 
averaged together to form a single image. This was done by strobing the function 
"Average Image Sequence” in the "Motion Sequence Capture/Display” sub-menu. The 
averaged image is stored in the buffer designated as the starting sequence buffer. 

The hard drive of the 80386 PC used in the present work contained software 
applications used in other research, which presented inadequate memory space for image- 
data storage. Therefore, all data was transferred from the frame-grabber board to floppy 
disk. The command to store data to disk 1s found in the main menu and is designated as 
"Image File Load/Save". From the resulting sub-menu, select the option "File Load/Save, 
Tiff Format w. AOI". The sub-menu for this option is brought up onto the computer 
screen. Position the cursor on the option, "Save Image to File. | Name", and type "b: 
filename.tif". Striking the "enter" key will bring up the menu "Select IMAGE BUFFER 
and/or AREA OF INTEREST" to appear; select the option "Image Area Of Interest: Full 
Image”. The image will then be sent to the floppy. "Full Image" refers to the image size 
that is set in the "Video Resolution" sub-menu; thus if the image resolution is set to 752 x 


480, this will be the image size that is stored on the disk. 
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Once ten averaged images have been saved to disk, they can be loaded back into ten 
frame buffers on the frame-grabber board using "Image File Load/Save" function. Again, 
select "File Load/Save, Tiff Format w. AOI" and then position the cursor on "Load Image 
from File. Name:". Type "b:filename.tif" and strike "enter". In the "Select IMAGE 
BUFFER and/or AREA OF INTEREST" sub-menu, select "Image Area Of Interest: Full 
Image" and strike "enter". The image will be loaded into the current image buffer. The 
current image buffer is selected by simply toggling the Pg Up/Pg Dn keys. 

With all ten averaged images loaded into the frame grabber, return to the main 
menu and select the "Motion Sequence Capture/Display" function. Insure the "Sequence 
Starting and Ending Buffers" reflect the sequence of buffers that the averaged images are 
stored in (See Figure C4). This complete, the images are averaged together using the 
"Average Image Sequence" option, the resultant image being stored in the buffer 
designated the sequence-starting buffer. This image represents either the wind-on and 
wind-off images that will eventually be ratioed or the averaged dark-current image that 
must first be subtracted from these two images. The three resultant images are again 
stored on floppy disk. 

2. Subtraction 

The averaged wind-on, wind-off, and dark-current images are loaded back into 
consecutive buffers of the frame-grabber board using the steps for loading files that were 
discussed above. From the main menu select "Image Processing", from which the image- 


processing sub-menu is displayed on the screen as in Figure C5. Choose the "Two Image 
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> Image Processing 


! COMMAND >> |« Quit Menu 
Histogram Displays |» Moments & Center of Mass 


TMS320 Resident Operations 
TMS320 Loadable Operations 


Simple Pixel Operations Morphological Operations 
Image Copy, Resize, Rotate Transform (FFT) Operations 
Two Image Arithmetic Image Interlace Modification 
Spatial Filters Equalization & Normalization 
Edge Detectors Two Image Normalizations 

NxN Operations 32 Bit Integration 
Miscellaneous Operations Image Sequence Operations 


Image Processing Modes 





Figure C5. "Image Processing" Sub-Menu 
Arithmetic" option, which sub-sequently displays the options available for conducting 
operations involving two images, this menu being illustrated in Figure C6. Four 
subtraction options may be selected; the option " Abs(PixB - PixA)" was used in the 
present work, where the absolute value of an intensity value is taken if it is less than zero. 
When this function is selected, the "Select IMAGE BUFFER and/or AREA OF 
INTEREST" for PixA menu is displayed. Before selecting anything from this sub-menu, 
ensure that the current-image buffer that is being displayed on the image monitor is that 
containing the dark current image; this represents PixA. From the "Select IMAGE 
BUFFER and/or AREA OF INTEREST" menu, choose the option "Full Area of Interest: 
Full Image". When this step has been completed, the "Select IMAGE BUFFER and/or 


AREA OF INTEREST" for PixB is displayed. Here, ensure the current image buffer 


122 


> Two Image Arithmetic 
I COMMAND >> |< Quit Menu 


! Add Images: PixB <- PixB + PixA MOD 256 
I Add Images: PixB Min(Pix + PixA, 255) 
I Subtract Images: PixB 128+(PixB-PixA)/2 

І Subtract Images: PixB PixB - PixA MOD 256 
I Subtract Images: PixB Max(PixB - PixA, 0) 
! Subtract Images: PixB Abs(PixB - PixA) 

! AND Images: PixB PixB & PixA 

! Exclusive OR Images: PixB PixB o PixA 

! Average Images: PixB (PixA + PixB) / 2 


I Product Images: PixB (c0*PixA+rc1) * (c2*PixB+c3) / có 
"Product Coefficient: cO 1 "Product Coefficient: c3 
“Product Coefficient: cl 0 “Product Coefficient: c4 
“Product Coefficient: c2 1 


I Ratio Images: PixB (с0:5.-Рдаанас д/д РЛХА 65) 
"Ratio Coefficient: cO 1 "Ratio Coefficient: c2 
“Ratio Coefficient: cl 0 “Ratio Coefficient: c3 


I Function on Images: PixB <- f(PixB, PixA). Expression: 


I UnInsert Images: If (PixB = PixA) PixB <- 0. Eps #: 
! Insert Images: If (PixB == 0) then PixB <- PixA 





Figure C6. "Two-Image Arıthmetic" Sub-Menu 
being displayed on the monitor is either the wind-on or wind-off image; this represents 
PixB. Again select "Full Area of Interest: Full Image" from the sub-menu. Execution will 
begin, the modified image residing in the frame buffer that originally contained PixB. 
Conduct this operation on both the wind-on and the wind-off images. 

3. Registration 

Registration of images will be unique to every application. As discussed in Chapter 

HI, the registration procedure between the wind-on and wind-off images for the jet/flat 
plate габой involved solid body motion, a relatively simple registration problem. 
Registration between images that involve more complex scenarios such as translational 
shifts, rotational shifts, and geometric distortions, will require more complicated software 


procedures than were used here. 
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If the extent of distortion between the images 1s simply solid body translation, then 
the following procedure can be used: Load the wind-on and wind-off images into 
consecutive buffers. From the main menu, select the "Pixel Peek, Poke & Plot" option; 
this sub-menu 1s shown in Figure C7. When this menu is displayed, the cursor will be 
overlaid on the image monitor, and can be controlled with the left button on the mouse is 
depressed. The cursor coordinates remain the same as the image buffer is changed, 
allowing examination of feature registration among image sequences. Thus using the 
wind-off image as the reference image, the vertical and horizontal characters of the cursor 
are aligned on a feature or control point that is distinctive in both images. By shifting back 
and forth between the buffers containing the images, using the Pg Up/Pg Dn keys and not 


disturbing the cursor location, an indication of the magnitude of distortion is obtained. 


» Pixel Peek, Poke & Plot 


! COMMAND »» |< Quit Menu 


! Cursor X Coordinate ! Current Image Buffer 
t Cursor Y Coordinate ! Device/Unit Selected 
! Pixel Value at Coordinate 


Peeking Within Image Buffers: Peeking Across Image Buffers: 
! Pixel Peek 17x17, X-Y I Pixel Peek 17x17, B-Y 

l! Pixel Plot, Line (X) ! Pixel Peek 17x17, X-B 

І Pixel Plot, Column (Y) ! Pixel Plot, 8uffers (B) 

|! Pi 





Figure C7. "Pixel Peek, Poke and Plot" Sub-Menu 
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Return to the main menu and select "Image Processing", followed by the option "Image 
Copy, Resize, Rotate"; this sub-menu is shown in Figure C8. The function "Copy & X,Y 
Shift" copies an image , shifting the image left or right and up or down in its buffer. The 
Left(-X) Right(+X) Shift parameter specifies the number of pixels the image 1s to be 
shifted left or night; the Up(-Y) Down(+ Y) Shift parameter specifies the number of pixels 
the image is to be shifted up or down. When the parameters have been entered, strobe the 
"Copy & X,Y Shift" function using the left button of the mouse or the "enter" key. 
(Ensure that the wind-on image is in the current image buffer prior to this step). The 
"Select IMAGE BUFFER and/or AREA OF INTEREST" for image A menu will appear 
from which the "Image Area of Interest: Full Image" Sption is selected; when this is done 


the "Select IMAGE BUFFER and/or AREA OF INTEREST" will re-appear for image B 


! COMMAND >> |< Quit Menu 


! Copy Image 


! Copy & X,Y Shift ! Copy & Flip Left/Right 
* Left(-X) Right(*X) Shift: ! Copy & Flip Top/Bottom 
" Up(-Y) Down(*Y) Shift: ! Copy & Flip Top/Bottom & Left/Right 


! Copy & Skew Left/Right ! Copy & Skew Up/Down 
T Skew at Top Edge: 0 ” Skew at Left Edge: 0 
7 Skew at Bottom Edge: 0 ^ Skew at Right Edge: 0 


Following operations require nonoverlapping source and destination: 


1 Copy & Resize м. Bilinear Interpolation 
! Copy & Resize w. Nearest Neighbor Interpolation 
7 Interpolate: Top->Up же) 
^ Interpolate: Top->Down ven) 
Interpolate: Top->Left гай! 
Interpolate: Top->Right С } 
” Interpolate: Тор Flipped NO 


! Copy & Resize м. Linear Area Interpolation 
! Copy & Resize w. Gaussian Area Interpolation 


! Copy & Resize w. Pixel Replication 
” Pixel Replication, X 
^ Pixel Replication, Y 2 


! Copy & Rotate м. Bilinear Interpolation. Angle: 
! Copy & Rotate w. Nearest Neighbor Interp. Angle: 
^ Rotation Correction for Aspect Ratio (X/Y): 





Figure C8. "Image, Copy, Resize, Rotate" Sub-Menu 
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where again the option "Image Area of Interest: Full Image" is chosen. When this 
selection is made the operation will then be executed. The current image buffer containing 
the wind-on image will be replaced with the modified wind-on image. Return to the "Pixel 
Peek, Poke & Plot" sub-menu and use the cursor to determine if the wind-on and wind-off 
images are properly registered. Several iterations may be necessary until the images are 
aligned. 
4. Ratioing 

With the images correctly registered, the ratioing process can now be executed, i.e., 
the wind-off may be divided by the wind-on images. Proceed to the "Two Image 
Arithmetic" sub-menu as previously done when the dark current image was subtracted 
from the wind-on and wind-off images. Under the "Ratio Images" option, set the ratio 
coefficients, that are used for "gain and offset", as follows: cO = 100 (or 128), cl = 0, c2 = 
1, с3 = 0. Next strobe the "Ratio Images" function which will cause the "Select IMAGE 
BUFFER and/or AREA OF INTEREST" menu for PixA to be displayed. Ensure the 
current image buffer is displaying the wind-on image. Select "Image Area of Interest: Full 
Image". This will cause the "Select IMAGE BUFFER and/or AREA OF INTEREST" to 
be re-displayed for PixB. Now shift the buffer containing the wind-off image such that it 
becomes the current image buffer. Again select "Image Area of Interest: Full Image". 
The ratioing process will then commence replacing the image designated as PixB, i.e., the 


wind-off image, with the modified image. 
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5. Thresholding 
A black background ıs usually added to the ratioed image to provide contrast and 
enhancement to the resultant image. From the main menu chose "Image Processing" and 
then the option "Simple Pixel Operations". Under the function "Threshold Pıxel Values" 


set the values as follows: 


- Threshold: Lowbound 0 
- Threshold: Highbound 255 
- Threshold: NewValue 0 


Strobe the "Threshold Pixel Values" function with the left mouse button; the "Select 
IMAGE BUFFER and/or AREA OF INTEREST" menu will appear. The boundary shape 
of the painted surface depicted in the ratioed image will dictate which area of interest to 
select for this image processing operation, 1.e., "Interactive Rectangular", "Interactive 
Elliptical", "Interactive Freehand", and "Interactive Polygon". The area of interest 
includes the area of the image that the painted surface does not occupy; for the flat plate, 
since the sides were parallel to the image axis, the "Interactive Rectangular Region" was 
selected as the region of interest. The area of the flat plate image not occupied by the 
plate was broken down to four areas of interest, that is, the area above the plate, below 
the plate, left of the plate, and to the right of the plate. Using the mouse with the right 
button depressed, these areas were outlined individually with the cursor. Once the areas 


had been outlined, the thresholding process was executed. 
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6. Pseudocoloring 
To add false coloring or pseudocoloring to the gray-scale image, select the 

"Contrast & Lookup Tables" option from the main menu. This menu is illustrated in 
Figure C9. Each of the primary colors, red, green, and blue are set by selecting the 
"Numerically Set & Show" option for the individual colors. The option displays the 
sub-menu illustrated in Figure C10. Using the function "Set Partial Linear Ramp" any 
configuration of linear ramps may be constructed. A graphical display of the inputs can be 
shown on the computer screen by selecting the "Plot Table" option in the menu illustrated 


in Figure C10. 


» Contrast & Lookup Tables 


|! COMMAND >> |< Quit Menu 


! Number of Lookup Table Sets 
|! Current Lookup Table Set 


Numerically Set & Show: Red Table 
Numerically Set & Show: Green Table 
Numerically Set & Show: Blue Table 


Interactive Contrast & Color: Center, Width, Depth 
Interactive Contrast & Color: Black Level, Contrast 
File Load/Save of RGB Set 





Figure C9. "Contrast & Lookup Table" Sub-Menu 
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> Define Lookup Table 


! COMMAND >> |< Quit Menu 


! Set Linear Ramp ! Set Partial Linear Ramp: 
! Set Log Curve “Starting Coordinate 

|! Set Exp Curve “Starting Value 

|! Set Gamma Curve, Gamma “Ending Coordinate 

1! Set as f(C). Expression: “Ending Value 


|! Complement/Invert Table 
1! Transpose Table 


? Show Table 
? Plot Table 


! Save Lookup Table to File. Name: 
! Load Lookup Table from File. Name: 





Figure C10. "Define Lookup Table" Sub-Menu 
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APPENDIX D. DARK CURRENT 


Thermal generation of electrons or "dark-current" is a naturally-occurring consequence 
of using semiconductor devices for photo reception. To understand dark-current effects 
on pressure-sensitive paint measurements, it is necessary to review the operation of a 
charged-coupled device (CCD). A schematic illustration of a CCD is shown in Figure D1. 

Consider a single CCD electrode in Figure D1. Absent the application of a bias to the 
gate electrode, a uniform distnbution of holes or majority carriers exists in the p-type 
semiconductor. As a positive voltage is applied to the gate, the holes are repelled from 
the semiconductor immediately beneath the gate, forming a depletion layer or potential 
well. An increase in the voltage causes the depletion region to extend further down into 
the bulk semiconductor, causing an increasing positive potential known as a surface 
potential to exist at the semiconductor and insulator (oxide) interface. Eventually a gate 
bias is reached at which the surface potential becomes so positive that the minority carriers 
(in this case electrons) are attracted to the surface where they form an extremely thin but 
very dense inversion layer. [Ref. 28:p. 6 ] 

In image-sensing applications of CCD's, the minority carriers are generated when 
photons entering the silicon substrate excite the valence electrons of the material into the 


“conduction band", becoming free electrons. The electrons are then collected under the 
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GATE OXIDE CCD ELECTRODES 


POTENTIAL WELL 
P-TYPE SILICON SUBSTRATE 


Desired event : weil filled with 
electrons as a result of photon energy 


Well "off-set" due to "Dark Current” 





Figure D1. Representation of the Effect of Dark-Current on the Creation of Charge 
electrode of the applied bias. The number of electrons under the electrode within a given 
period of time is proportional to the local light intensity. 

In addition to the photoelectric effect described above, free Р can also be 
produced by thermal energy, of which the rate of generation increases exponentially with 
temperature for a pure semiconductor. Thus, the presence of thermally generated minority 
carriers limits the maximum possible usefulness of the potential well by displacing optically 


generated carriers. This produces an "offset" in the measurement image such that the 
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zero-incident light level does not correspond to a zero gray-level pixel index. Thus, the 
dark image must be removed from the measurement images. To do this, a dark image is 
acquired by prohibiting any light from striking the CCD array and sampling an image. 
Several of these images are taken and time averaged together. This is necessary because 
dark current generation tends not to be uniform over a whole array. The resulting 
averaged dark current is then eliminated from the measurement images using image 


processing techniques discussed in Chapter III. 
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APPENDIX E. SUPERSONIC NOZZLE DESIGN PROGRAM 


Software for the design of two-dimensional supersonic nozzles was acquired from the 
Naval Surface Warfare Center (NSWC), White Oak, Maryland. The software package is 
comprised of two programs; the first program is concerned with establishing the 
isentropic contour of the nozzle while the second corrects for the viscous effects of the 
fluid and adjusts the contour accordingly. 

Support documentation for the software from a user's standpoint is not available. In 
addition, the early fortran code is not well structured, making it difficult to interpret the 
algorithm and follow its logic. It is difficult to recognize many of the vanables, and there 
are no comment statements. However, enough information has been deduced such that 
the program may be operated effectively. The intention of this appendix is to document 
what is now known about the software. 

Modifications have been made to the original software primarily involving the manner 
in which program input entries are made and data storage accomplished. Originally, the 
program was designed to access punch-card readers and magnetic-tape drives for data 
input and output management. Entries are now made directly to the program through 
interactive software in a format that is more flexible and user friendly. Output data are 
stored in a file structure compatible with hard drive systems. A listing of the modified 


software is given in Section C of Appendix E. 
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A. INVISCID CONTOUR SOFTWARE (Program char6) 

The program "char6" generates the isentropic contour of the supersonic nozzle using 
the method of characteristics. Initiation of the program will depend on the filename 
assigned to the object code at the time of compiling and linking of the source code. Once 
started, the program first displays an introductory summary as shown in Figure El. As 
indicated by these remarks, the storage of the output data is currently limited to the data 
for one run. A successive run will overwrite the data storage files from the previous run. 

Entries made by the user through the interactive sequence of prompts for input data 
have been sub-divided into two categories; the first category identifies the "characteristic" 
lengths and nozzle type that will be used to define the configuration; the second specifies 


dimensions and flow characteristics to which the nozzle will be designed. 


U.S. NAVAL ORDANCE LABORATORY, WHITE OAK, MARYLAND 
COMPUTER PROGRAM FOR THE DESIGN OF TWO-DIMENSIONAL SUPERSONIC NOZZLES 
PART | 


THE ISENTROPIC CORE OF THE NOZZLE 


INTRODUCTION: 


This program uses the method of characteristics to determine the contour of a two-dimensional 
isentropic core of a supersonic nozzle. The program requires two sets of inputs; one that 
provides the general characteristic lengths and nozzle types that will be used to describe the 
nozzle; the second set outlines the specific flow characteristics and geometric lengths that the 
contour will be designed to. The program is currently limited to one output file; in other words 
successive runs will overwrite data generated on previous runs. 





Figure El. Illustrative Example of the Introductory Remarks for Program char6 
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1. Category 1: Characteristic Lengths and Nozzle Types 


1. INPUT A DESCRIPTION OF THE PROJECT 


Allows a short summary of the project to be recorded as part of the final output 
summary. The description ıs limited to 80 characters. 


2. INDICATE TYPE OF NOZZLE, 1-SHORT, 2-REGULAR, 3-EXPONENTIAL 


Provides the option of choosing a specific contour shape. The methodology the 
program uses to generate the contour begins by first determining the centerline 
Mach number distribution based on the choice of contour. The method of 
characteristics is then used in conjunction with the centerline distribution in 
determining contour geometry. 


3. INDICATE THROAT OR EXIT HEIGHT TO BE READ, 1-ЕХІТ, 
2-THROAT 


Irrespective which is used; the program calculates the option not selected using 
simple isentropic relations. Both heights are employed in the program in 
computing the contour. 


4. INDICATE CHARACTERISTIC HORIZONTAL LENGTH TO BE READ, 
1-EXIT LENGTH, 2-START OF TEST RHOMBUS 


Refer to Figure E2 for clarification of the honzontal lengths discussed. The exit 
length is defined as the distance between the nozzle throat and the nozzle exit. 


The characteristic network created by the method of characteristics 1s such that 
the last two characteristic lines in the grid complete the expansion process by 
turning the flow parallel to the centerline of the nozzle and thus expanding it to 
the desired test section Mach number. These two characteristic lines can also be 
thought of as forming half of an imaginary rhombus as shown in Figure E2. This 
is known as the test section rhombus. The distance between the nozzle throat 
and the intersection formed by the two characteristic lines defines the horizontal 
length referred to above. 


5. INDICATE OUTPUT, 1-STEP-BY-STEP, 2-WALL POINTS ONLY 


Generally, the approach used in applying the method of characteristics can be 
summarized in a three step process: Step 1 -- the determination of the 
characteristic lines in the xy space using 
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Test Section Rhombus \ 
\ 


Sonic Line 


Throat Height 


(v2 
1 
Initial Data Line 
Start of Test Rhombus == 


| Exit Length 
Figure E2. Schematic of Supersonic Nozzle Design by the Method of Characteristics 





(Donar =tan(0 + 11) 


where 9 is the angle between the streamline of the flow and the horizontal and 
и 15 the angle between the streamline and the characteristic line; Step 2 -- the 
determination of the compatibility equations 


O+v(M) =const=K_ (along the right characteristic) 
0 — v(M) = const = К+ (along the left characteristic) 


where again 9 is the angle between the horizontal and the streamline, and, v(M) 
is the Prandtl-Meyer function, the compatibility equations holding along the 
characteristics; Step 3 -- the solution of the compatibility equations point by 
point along the characteristics. 


Note: The initial data line used in starting the computation 
of the method of characteristics in this program begins with 
the assumption that the sonic line at the throat is straight. A 
point AM downstream of the sonic line along the centerline 
Mach number distribution is chosen as the initial point for the 
characteristics procedure. With the combination of this 
initial point , the centerline distribution, and the radius of 
curvature at the throat, the method of characteristics 1s carried 
out in a marching fashion. 
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Choosing the option step-by-step will instruct the program to include in the final 
output all "steps" used in calculating the characteristic net, including local data 
pertaining to the nodal points comprising the turning angle, the x and y coordinates 
of the grid points and the Mach number. If in turn the option "wall points only” is 
exercised, the data associated with the grid points that decide the contour or wall 
geometry are included in the final output. 


Once the last selection for the first category of inputs is chosen and the "enter" key is 
subsequently struck, the selections made are then displayed in the following format: 
NOZZLE = SHORT (REGULAR, EXPONENTIAL) 
HEIGHT TO BE READ = EXIT (THROAT) 
LENGTH TO BE READ = LENGTH TO EXIT (LENGTH TO RHOMBUS) 
OUTPUT FORMAT = STEP-BY-STEP (WALL POINTS) 
The question, "Are input values correct?, 1-Yes, 2-No", follows immediately, allowing the 
opportunity to make changes to the inputs if necessary. Selecting "Yes" causes the 
following menu to be displayed: 
INPUT VALUES TO BE CHANGED (1, 2, 3, 4): 
1- ТҮРЕ OF NOZZLE: A. SHORT, B. REGULAR, C. EXPONENTIAL 
2- HEIGHT TO BE READ: A. EXIT, B. THROAT 
3 - LENGTH TO BE READ: A. EXIT, B. TEST RHOMBUS 
4-TYPE OF OUTPUT: A. STEP-BY-STEP, B. WALL POINTS ONLY 
By selecting a value from one to four a sub-menu of the descriptors to the nght of the 
colon of the input chosen is listed. For example, selecting the value of 2 (Height to be 
read) from the above menu yields 
1 - EXIT, 2- THROAT 
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prompting a choice of height designation. After hitting the "enter" key, the program will 
again display the first-category input selections, showing the latest modifications. An 

opportunity to make additional alterations is also afforded when the question, "Any other 
changes?, 1-Yes, 2-No", immediately follows. Selection of "No" initiates the sequence of 


events for data input associated with the second category. 


2. CATEGORY 2: NOZZLE FLOW AND GEOMETRY SPECIFICATIONS 
Commencement of the sequence of prompts for the second category of inputs is 


denoted by the statements 


"THE NEXT SET OF INPUTS SPECIFY FLOW CHARACTERISTICS AND 
SPECIFIC GEOMETRY DIMENSIONS FOR THE NOZZLE" 


***** PLEASE UTILIZE DECIMAL POINTS FOR THESE INPUTS ***** 


It is necessary to use decimal points when specifying numeric data; the program is 
formatted to read the input data in scientific notation, though numeric values may be 
entered into the program as real numbers. If a decimal point is not explicitly stated, the 
exponent field in the scientific notation format will reflect this with zeros, resulting in an 
ambiguous input. Thus, if the real number "14" is the desired input, it must be entered as 
"14."; this results in the scientific notation representation "0.14000E+02". Ifthe number 
is entered as "14", the resulting representation will be "0.14000E+00". 

The inputs for the second group are outlined as follows: 

1. INPUT RATIO OF SPECIFIC HEATS 


Self-explanatory. 
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2. INPUT TEST SECTION MACH NUMBER 


The Mach number of the uniform flow that is required at the exit plane of the 
nozzle. 


3. INPUT EXIT LENGTH (INCHES) OR INPUT BEGINNING OF TEST 
RHOMBUS (INCHES) 


The choice made in the first category of inputs concerning the horizontal length will 
dictate which prompt will be displayed. Refer to Figure E2. 


4. INPUT HALF HEIGHT OF EXIT (INCHES) OR INPUT HALF HEIGHT OF 
THROAT (INCHES) 


Again dependent on the selection made in the first category of inputs with 
reference to the height to be read. 


5. INPUT RADIUS OF CURVATURE AT THROAT (INCHES) 
Self-explanatory. 
6. INPUT MULTIPLICATIVE FACTOR FOR "X" INCREMENT 


Influences the spacing between consecutive contour points. A range between 1.0 
(representing 150 points to describe the contour) to 0.5 (representing 400 points) is 
available. 


Note: Generation of the cubic centerline Mach number distribution, 
which is associated with the regular nozzle contour, is dependent 
upon the result of a "feasibility expression" involving the second 
category values input by the user. If the outcome of the expression 
does not fall into a pre-determined interval, the program halts 
execution and displays the message, 


*** CUBIC CENTERLINE MACH NO. DISTRIBUTION IS 
NOT APPROPRIATE FOR INPUT GIVEN *** 


At this point, the program must be re-initialized and the process of 
re-entering the input started over. One of the input values, gamma, 
test section Mach number, exit length or length to test rhombus, exit 
or throat height, or the radius of curvature, or a combination thereof 


14] 


must be changed such that the compatibility requirements set forth by 
the feasibility expression are satisfied. 


When the data for the multiplicative factor has been entered and the "enter" key struck the 


second category of inputs are displayed in the following fashion: 


GAMMA = 0.14000E+01 
SMT = 0.17000E+01 
XEXXT = 0.11500E+02 
YTHEX = 0.19500E+01 
RCT = 0.11500E+01 
DELMUL = 0.70000E+00 
ARE INPUT VALUES CORRECT ?, 1 - YES, 2 - NO 


where 


GAMMA = Ratio of specific heats 
SMT = Test section Mach number 

XEXXT = Exit length or length to the beginning of test 
rhombus depending on the choice in first category 
of inputs 

YTHEX = Half height of exit or half height of throat depending 
on the choice in the first category of inputs 

DELMUL = Multiplicative factor for "x" increment 


As in the first category of inputs, the option to make changes to the input data 15 


available. When changes are desired a listing of the input descriptors is displayed: 


INPUT VALUE TO BE CHANGED: 1 - GAMMA, 2 - SMT, 
3 - XEXXT, 4 - YTHEX, 5-RCT, 6 - РЕМИ, 


A selection will exhibit a directive prompting the user to input the new value. For 
example, choosing gamma (1) results in, 

INPUT NEW VALUE FOR GAMMA 
Fulfilling this and subsequently striking the "enter" key causes the second category of 
inputs to be re-displayed showing the recent modifications. Additional changes may also 


be conducted at this point if necessary. 
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If the second category of inputs have been reviewed with satisfaction, the program will 
begin calculation of the isentropic contour at this point. Conclusion of the routine is 


indicated by the remarks, 


OUTPUT ISENTROPIC DATA FOR NOZZLE CONTOUR 
STORED IN FILENAME "outpt1" 


DATA FROM THIS FILE TO BE USED IN THE VISCOUS 
PROGRAM "пЫб.Ғ" IS STORED IN FILENAME "outpt2" 


The file "outpt1" contains a systematic breakdown of the data calculated by the program 
which are then available to the user for review. A complete explanation and an example of 
the format are presented in the following section. The file "outpt2" contains the data that 
are to be utilized in the program "nbl6" and are in a format that is convenient for file 
processing. 

3. "outptl" EXAMPLE AND DESCRIPTION 

Figure E3 shows the format in an example of the file "outpt1". It has been 

sub-divided into it's major components and labeled for reference below. 


1. Introductory Header - Self-explanatory. 


2. Descriptive Comments - Refers to the comments made in the first category 
of inputs when prompted by the program for a description of the project. 


3. User Inputs - Encompasses all the inputs into the program made by the user. 
The input identifiers associated with the first category of inputs are defined 
as follows: 


INDSR = Nozzle Type 

INDYTH - Characteristic Height (Throat or Exit) 
INDXTE - Characteristic Length (Exit or Test Rhombus) 
INDOUT - Output Mode (Step-by-Step or Wall Points) 
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U.S. NAVAL ORDANCE LABORATORY, WHITE OAK, MARYLAND 


COMPUTER PROGRAM FOR THE DESIGN OF TWO-DIMENSIONAL SUPERSONIC NOZZLES INTRODUCTORY 


HEAD 
PART I ER 


THEISENTROPIC CORE OF THE NOZZLE 


1.7 MACH NOZZLE DESCRIPTIVE 
COMMENTS 


INPUT DATA 
CARD NO. 1 INDSR = 2 
INDYTH = 1 
INDXTE = 1 
INDOUT = 1 
INDCAL= 0 
INDTAP= 0 
INDCRD= 0 USER 
INPUTS 
CARD NO. 2 GAMMA = 0.14000E+01 
SMT = 0.17000E+01 
XEXXT = 0.11500E+02 
YTHEX = 0.19800E+01 
RCT = 0.11500E+02 
DELMUL= 0.70000E+00 
NOZZLE PARAMETERS 
MACH NO. YTH YEXIT XEXIT XT NOZZUE 
0.17000E+01 0.14803E+01  0.19800e+01 0.11500E+02 0.87780E+01 PARAMETERS 
RCT B 
0.11500E+02 0.35425E+00 
CENT ERLINE MACH NUMBER DISTRIBUTION USING CUBIC EQUATION 
N x MACH NO. 
1 0.534600E-03 0.100010E+01 CENTERLINE 
2 0.732600E-03 0.100014E+01 DISTRIBUTION 
3 0.930600Б-03 0.100017Е+01 
CHARACTERISTICS MANIPULATION (‘INTERPOLATED WALL VALUES) 
N K X(INCHES)  Y(INCHES) MACHNO. STREAM ANGLE 
2 3 0.629684E-03 0.6428261E-02  0000000Е-00 0.136713E-04 MET HOD OF 
3 2 0.828632Е-03 05607560Е-02 0.000000E+00 0.156680Е-04 CHARACTERISTIC 
3 3 0.718964E-03 0.1198140E-01 0.000000Ё-00  0.293399Б-04 RESULTS 
FINAL OUTPUT TO PART |, ISENTROPIC CORE CONTOUR 
| ХЇМСНЕ5) || У(ІМСНЕ5) CLMACHNO. WALL MACHNO. WALLANGLE RAD OF CURV (IN) CONTOUR 
1 0.435619E+00 0.149273E*01 0.107887E+01 0.114198E+01 0.220395E+01 0000000Е-00 OUTPUT 


2 0.462758Е-00 0.149380E+01 0.108359E+01 0.114598E*01 0.228777E+01 0. 188330Е-02 
З 0.490391Е+00 0.149492Е+01 0.108838E+01 0.115003E*01 0.237073Е+01 0.193154Е+02 


Figure E3. Example Showing the Format of File "outpt1" 
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INDCAL, INDTAP, and INDCRD are no longer user-accessible. These 
variables were used in the original software and controlled the medium to be 
used to record the output data, i.e., magnetic tape or punch cards. The 
variables have been "hardwired" in order to accommodate the current 
software modifications for data output management. 


The input descriptors for the second category of inputs have been previously 
defined. 


4. Nozzle Parameters - These include, in addition to those dimensions 
supplied by the user, the "related" values calculated by the program using 
isentropic relationships or simple geometric relations. Thus, given the input 
of the throat height the program calculates the exit height ,or, given the 
horizontal length from the throat to the exit plane the program calculates the 
length from the throat to the test rhombus. The parameters are defined as, 


Mach No. = Nozzle exit Mach number 
YTH = Half height of nozzle throat 
YEXIT = Half height of nozzle exit 
XEXIT = Length between nozzle throat and exit 
XT = Length between nozzle throat and test rhombus 
RCT = Radius of curvature at throat 
* B = Indicates result of "Feasibility Expression" 


* Note: The result of the feasibility expression, to determine 
if the cubic centerline Mach number distribution may 

be calculated with the given nozzle dimensional inputs, 

must fall in the interval -2.0 « B « 1.0. 
5. Centerline Mach Number Distribution - Denotes the centerline location 
(relative to the throat) and corresponding Mach number. The function used 
to generate the distribution, i.e., cubic, exponential, etc., is dependent upon 
the choice of nozzle type (regular, exponential, short). 


6. Method of Characteristic Results - Describes the characteristic manipulation 
used to generate the wall values along the contour. 


N = Characteristic line (From left to right) 
L = Nodal points along characteristic line 


X (inches) = X-coordinate of nodal point 


145 


Y (inches) = Y-coordinate of nodal point 
Mach No. = Local Mach number at nodal point 


Stream angle = Local velocity angle at nodal point with 
respect to the horizontal 


7. Contour Output - Final output of the isentropic core. 
X(INCHES) = X - coordinate of isentropic core 
Y(INCHES) = Y - coordinate of isentropic core 
CL MACH NO. = Centerline Mach number 


WALL MACH NO. = Mach number along contoured wall 
WALL ANGLE = Wall angle relative to centerline 


RAD. OF CURV. (IN) = Radius of curvature 


B. MODIFICATION FOR VISCOUS EFFECTS (Program nbl6) 


The program “nbl6" modifies the contour of the nozzle to account for the boundary 
layer formation. Originally, this software considered the boundary layer effects along the 
contour of the nozzle only and neglected the effects that would be caused with the 
addition of side walls. To compensate for this deficiency in correcting for boundary layer 
displacement, the approach used by Demo [Ref. 29: p. 24] was employed. The 
assumption is made that the displacement thickness is the same on both contoured and 
side walls. The viscous displacement effects on all surfaces are then summed together and 


evenly distributed between the two contour surfaces only, as illustrated in Figure E4. 
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Figure E4. Nozzle Wall Correction for Boundary Layer 


Thus, the total displacement area of the entire nozzle cross-section is compensated by 
modifying only the contoured walls, allowing the side walls to remain flat and parallel. 
The displacement required to the inviscid contour was calculated using the expression 
S (El) 

where 0” = displacement thickness, B = nozzle width, and y = half height of the contour. 
This equation has been incorporated into the software. The final output, discussed later, 
was modified to reflect this change. 

The data-output management structure, as before, is limited to storing data for one 
run. A successive run will overwrite the data stored from the preceding run. 

Input into the program is supplied by both the output from "charó" (1.e., file outpt2) 


and the user using an interactive presentation similar to that discussed above for the 
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inviscid contour. The first category of inputs consists of a series of functional mode 
selections concerned with the thermal and velocity boundary layer behavior. The second 
category is comprised of numerical input prompts for flow conditions and nozzle 


dimensions. 
1. CATEGORY 1: FUNCTIONAL MODES 
The functional modes establish the desired profile behavior that the velocity and 
thermal boundary layers will undergo while traversing the length of the nozzle. The 


modes are outlined below. 


1. INDICATE WHERE THE BOUNDARY LAYER CALCULATION 
WILL TAKE PLACE, 1. ALONG CENTERLINE ONLY, 2. ALONG 
THE WALL ONLY, 3. BOTH THE CENTERLINE AND THE WALL 


Specifying the first option results in the message, 


Ж Ж Ж Ж Ж Ж ж ж ж Ж ж Ж ж Ж жж Ж e e e oe e e oe e e ЖҰЖ e e e e o ЖЖ 


" CENTERLINE PARAMETER INPUT DATA T 

жжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжж ж 
indicating that all input into the program from this point will pertain to 
the effects of the boundary layer along the nozzle centerline. Selection 
of the second option displays 


Ж Ж Ж ж Ж Ж Ж ж ж Ж ж ж ж Ж Ж Ж ж ж Ж ж ж ж ж Ж Ж Ж ж Ж ж ж ж ж ж ж ж 


А WALL PARAMETER INPUT DATA Ё 


Ж Ж Ж Ж Ж Ж ж ж ж Ж ж Ж ж ж ж ж ж ж Ж Ж ж Ж ж Ж Ж ж Ж ж ж Ж ж ж ЖЖ Ж 


indicating further inputs will apply to the effects of the boundary layer 
along the nozzle contoured walls. Option three effectively causes 

the program to run twice, once for the effects along the centerline 
followed by a run associated with the boundary layer effects along the 
contoured walls. The final output for option three is consolidated as 
a single output vice individual outputs for each run. 
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2. INDICATE BOUNDARY LAYER MODE BEGINNING AT NOZZLE 
THROAT, 1. STRICTLY LAMINAR, 2. STRICTLY TURBULENT, 
3. LAMINAR-TRANSITION-TURBULENT 


Self-explanatory. 


3. INDICATE WALL TEMPERATURE FOR NOZZLE, 1. CONSTANT 
TEMPERATURE, 2. ADIABATIC WALL, 3. TEMPERATURE 
DISTRIBUTION ALONG NOZZLE TO BE INPUT BY PROGRAM 
USER 


Wall temperature is assumed only along the contour of the nozzle. Since 
the boundary layer displacement thickness for the entire cross-section of 

the nozzle is a function of the displacement thickness generated along the 
contour, it can be thought of that the temperature distribution exists both 
on the contour and side walls. 


4. INDICATE THE CHARACTERISTIC PARAMETER ON WHICH 
TRANSITION WILL BE BASED, 1. REYNOLDS NUMBER, 2. "X" 
DISTANCE FROM THE THROAT 


The transition mode only becomes active if the laminar-transition- 
turbulent selection of the boundary layer mode is chosen. The point 
of transition from laminar flow to turbulent flow may be based upon a 
horizontal distance "x" relative to the nozzle throat and along the 
centerline or a Reynolds number in which the momentum thickness of 
the boundary layer is used as the characteristic length 


The last entry made for the mode options is followed by a listing of the selections to 
verify the settings: 


B.L. MODE = LAMINAR (TURBULENT, LAM-TRANS-TURB) 


WALL TEMP. MODE = CONSTANT (ADIABATIC, USER 
DEFINED) 


*TRANSITION PARAMETER = REYNOLDS NUMBER ("X" 
DISTANCE) 


B.L. CALCULATION = CENTERLINE ONLY (WALL ONLY, 
BOTH CENTERLINE AND WALL) 
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"ARE INPUT VALUES CORRECT”, 1-YES, 2-NO 
*Note: Displayed only if B.L. Mode = Lam-Trans-Turb 
If a mode is incorrectly set or a change is desired, "no" is subsequently chosen, resulting in 
a modification menu: 
INPUT VALUE TO BE CHANGED (1,2,3,4): 
l. B.L. MODE: A. LAMINAR, B. TURBULENT, C. LAM-TRANS-TURB 


2. WALL TEMP. MODE: A. CONSTANT, B. ADIABATIC, C. USER DEFINED 


3. B.L. CALCULATION: A. CENTERLINE, B. WALL, C. CENTERLINE 
AND WALL 


*4. TRANSITION: A. REYNOLDS NUMBER, B. "X" TRANSITION 
*Note: selection of the transition mode with the B.L. Mode 
set other than B.L. Mode = Lam-Trans-Turb will result in 


the message, "Transition parameter not required for inputs 
given". 


The mode to be modified is selected from the menu list, which in turn causes a sub-listing 
of the options to which the mode may be set to be displayed. For example, if the 
boundary layer mode is to be changed, its selection from the modification menu would be 
followed by, 

1. LAMINAR, 2. TURBULENT, 3. LAM-TRANS-TURB 
A choice of one of the available options sets the mode and then causes a re-display of all 
the modes and their current status; 

B.L. MODE = LAMINAR 


WALL TEMP. MODE = CONSTANT 
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TRANSITION PARAMETER = REYNOLDS NUMBER 
B.L. CALCULATION = WALL ONLY 
"ANY OTHER CHANGES?, 1-YES, 2-NO" 
Other modifications may be done at this point in a sequence similar to that discussed 
above. If no other changes are required or if the modes initially set were correct, the 
series of prompts for the second category of inputs will be initiated. 
2. CATEGORY 2: FLOW CONDITIONS AND NOZZLE DIMENSIONS 
The second category of inputs includes the numerical data associated with the flow 
conditions, and additional dimensions of the nozzle, which will be made clear below. The 
onset of the second set of inputs is denoted by the message 
* * * PLEASE UTILIZE DECIMAL POINTS FOR THE FOLLOWING INPUTS * * * 
followed by the first prompt for information. Decimal points are necessary for the 
numeric inputs for the reasoning previously explained. 
Most of the input prompts are self-explanatory and require no elaboration. Those 


prompts requiring explanation or having relevant comments will be summarized following 


the prompt. 


1. INPUT REYNOLDS NUMBER FOR TRANSITION 
(Reynolds number based on momentum thickness of boundary layer) 


Displayed only if the boundary layer mode is set to B.L. Mode = 


Lam-Trans-Turb and the transition indication mode set to Transition 
Parameter = Reynolds Number. 
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10. 


11. 


12, 


. INPUT TRANSITION DISTANCE FROM THE THROAT IN FEET 


Displayed only if the boundary layer mode is set to B.L. Mode = 
Lam-Trans-Turb and the transition indication mode js set to Transition 
Parameter = "X" Distance From Throat. 


. INPUT MOMENTUM THICKNESS OF BOUNDARY LAYER AT THE 


THROAT IN FEET 


A value of 1.0E-5 feet is suggested. 


. INPUT STAGNATION PRESSURE IN LBS/IN’ 
. INPUT TOTAL TEMPERATURE OF SUPPLY MEDIUM IN DEGREES R 


. INPUT TEMPERATURE OF THE WALL IN DEGREES R 


Displayed only if the wall temperature indication mode 15 set to Wall 
Temp. Mode = Constant. 


INPUT THE EXTENSION LENGTH OF NOZZLE PAST THE EXIT IN 
INCHES AS DEFINED FOR THE ISENTROPIC CORE 


The nozzle contour may be extended past the exit plane of the nozzle, 
as defined by "char6", if, for instance, a test section is to be added. 
Assumes uniform duct flow. 


. INPUT INCREMENT OF "X" IN INCHES FOR EXTENSION LENGTH TO 


INTERPOLATE CONTOUR OF NOZZLE 


Establishes the horizontal distance between consecutive points that 
describe the contour of the extension length. 


. INPUT WIDTH OF NOZZLE 


INPUT RATIO OF SPECIFIC HEATS 
INPUT EXPONENT FOR VISCOSITY-TEMPERATURE RELATION 
For air, this is (O —) 0.76. 


INPUT MOLECULAR WEIGHT OF GAS SUPPLY (AIR = 23.97) 
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13. INPUT SPECIFIC HEAT OF GAS IN BTU/LB-DEG-R (AIR = 0.24) 


14. INPUT PRANDTL NUMBER FOR GAS (AIR = 0.72) 


Following the conclusion of the last prompt, a listing of all the input descriptors and the 


values assigned to them is displayed: 


*REYNOLDS NUMBER = 0.50000E-02 

*"X" TRANSITION = 0.60000E+01 
MOMENTUM THICKNESS = 0.10000E-04 
STAGNATION PRESSURE = 0.47500E+02 
STAGNATION TEMPERATURE = 0.68000E+02 

**WALL TEMPERATURE = 0.72000E+02 

EXTENSION LENGTH PAST NOZZLE EXIT = 0.12000E+02 
INCREMENT "X" FOR EXTENSION = 0.65000E+00 
WIDTH OF NOZZLE = 0.40000E+01 
GAMMA = 0.14000E+01 
VISCOSITY-TEMP. EXPONENT = 0.76000E+00 
MOLECULAR WEIGHT OF GAS = 0.28970E+02 
SPECIFIC HEAT = 0.24000E+00 
PRANDTL NUMBER = 0.72000E+00 


"ARE INPUTS CORRECT?, 1-YES, 2-NO" 


*Note: Displayed only if boundary layer mode is set to B.L. 
Mode = 3 and the transition indication mode is set to 
either Transition Parameter = Reynolds Number or 
Transition Parameter = "X" Distance from Throat 


**Note: Displayed only it the wall temperature indication 
mode is set to Wall Temp. Mode = Constant 


If the option is exercised to change any of the input values, the user will be presented with 


the following menu: 


INPUT VALUE TO BE CHANGED (1,2,3,4,5,6,7,8,9, 10, 11, 12,13): 


. TRANSITION PARAMETER 
MOMENTUM THICKNESS AT THROAT 
. STAGNATION PRESSURE 

. TOTAL TEMPERATURE 

. WALL TEMPERATURE 
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л & YN = 


6. EXTENSION LENGTH OF NOZZLE 
7. "X" INCREMENT FOR EXTENSION LENGTH 
8. WIDTH OF NOZZLE 
9. RATIO OF SPECIFIC HEATS 
10. VISCOSITY-TEMPERATURE EXPONENT 
11. MOLECULAR WEIGHT OF GAS 
12. SPECIFIC HEAT (Cp) FOR GAS 
13. PRANDTL NUMBER 


Selection from the menu generates a prompt for the new information to be input for the 
chosen descriptor. Upon entering the new value and striking the "enter" or "return" key, 
the list containing the input descriptors and their assigned values will be re-displayed 
exhibiting the latest changes and soliciting further modifications. 

If at this point no other changes are desired and if the wall temperature indication mode 
is set to an option other than Wall Temp. Mode = Temperature Distribution Input By 
User, the program will proceed to calculate the modified contour according to the inputs 
given. If the user has elected to enter a defined temperature distribution along the nozzle 


contour, the following sequence of prompts can be expected: 


INDICATE THE NUMBER OF TEMPERATURE POINTS 

TO BE ENTERED. THREE OR MORE POINTS ARE 
REQUIRED OR PROGRAM WILL TERMINATE THE DATA 
OUTPUT PRIOR TO COMPLETION. 


INPUT THE "X" LOCATION FOR POINT NUMBER | 
RELATIVE TO THE THROAT 


INPUT TEMPERATURE IN DEGREES R FOR POINT 
NUMBER _ 


The second and third statement are repeated until the number of points designated in the 


first statement have been entered. When the last point is entered, the program proceeds in 
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its calculation of the modified contour. Completion of the routine is denoted by the 


message, 


PROGRAM RUN COMPLETE. OUTPUT DATA 
TRANSFERRED TO FILE "outpt3" 


3. FILE "outpt3" FORMAT DESCRIPTION 
Figure E5 illustrates an example of the format used for the output file "outpt3". 
Again it, as before, has been sub-divided into it's major components and labeled as 
follows: 
]. Functional Modes 
*INDINP = 2 


*INDPRG = 1 


INDMO=2 (Boundary Layer Mode) where 1 = Laminar, 
2 = Turbulent, 3 = Lam-Trans-Turb 


INDTW= 1 (Wall Temperature Mode) where 1 = Constant Temp., 
2 = Adiabatic, 3 = User Defined Distribution 


INDTR =1 (Transition Indication Mode) where 1 = Transition at 
Reynolds Number, 2 = Transition at "X" Length 


*INDCF = 1 
*INDOUT = 0 


*Note: Indicates values are hardwired into the program and 
are no longer user accessible 
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INPUT DATA ~ 


INDINP = 2 
INDPRG= 1 
INDMO = 2 
CARD 1 INDTW = 1 FUNCTIONAL 
INDTR = 0 MODES 
INDWC = 3 
INDCF = 1 
INDOUT = 0 
TRCON = 0.000000€ +00 
THZERO = 0. 100000€-04 
РО з 0.728000£+02 
CARD 2 TO = 0.820000E +03 
TWALL = 0.530000E +03 
XTEND = 0.285000E+02 OW CONDITIONS 
DL - 0.825000E-01 AND NOZZLE 
DIMENSIONS 
GAMMA = 0. 140000€+01 
OMEGA = 0.760000E+00 
CARD 3 WATE = 0.290000Е+02 
CP = 0.240000Е-00 
PR = 0.7 20000E+00 
XINC = 0.000000E+00 
* WALL DATA * 
X(IN) Y(IN) M 
0.318238Е-00 0.142340E +01 0.111979E+01 INPUT DATA 
0.333082E+00 0.142386E+01 0.112352Е-01 
0.348089E+00 0.142434E+01 0.112726E+01 
* OUTPUT CODE *- 

X(IN) S(IN) Y (IN) M REL(/FT) RETH PE/PO 
UE(FT/SEC) MUE(PSEC/FT2) RHOE(PSEC2/FT4) TER) TMR) TAD(R) ТЕЛО 
DELTA(IN) THETA(IN) DELSTAR(IN) H тато TWTAD TAD/TO 

DPOS(PIFT2) OMDS(1/FT) DUDS(1/SEC) DTHDS ТАШАЦР/ЕТ2) СЕ ВЕТАТН 
Q(BTU/FTZSEC) HA(BTU/FT2SECR) N 
OUTPUT 
POINT NO. 1 CODE 
0.263755E+02 0.283087Е402 0.1896006+01 0.170000E+01 0.120507E+08 0.300775E+05 0.202563E+00 
0.189944E+04 0.374381Е-06 0.237520E-02  0.519645E+03 0.530000€+03 0.788847E_03 0633714Е400 
0.348474E+00 0.300135E-01 0540098Е-01 0.160251E+01 0.646341Е%00 0.671866E+00 0.962008E+00 
0.000000E+00  0.000000€+00  0.000000£+00 0.100080£-02 0.857629Е4-01 0200005Е-02 0.000000E+00 
0.112397Е402 0.434220Ё-01  0.753040E «01 
0.307368E+02 
POINT NO. 2 
0.284380Е%02 0.284612E+02 O189600E+01 0.170000E+01 0.120507E+08 0.301403Е+05 0.202563Е-00 
0.189944Е-04 0.3743815-06 0.237520E-02 0519645£+03 0.530000£+03 0.788847E+03 0.633714E+00 
0.349255E+00 0.200761Е-01 0.542106E-01 0.180245Е-01 0.646341Е-00 0671886Е-00 0.98200BE+00 
0.000000E+00 0.0000006+00 0.000000=+00 0100047E-02 0.857346E+01 0200095Е-02 0.000000E+00 
0.112397Е+02 0.434220E-01 0.753225Е+00 
0.307953E «02 
* WALL RESULTS * 

X(INCHES) DELSTAR(INCHES) YCORE+DELSTAR MWALL DELTAY YOORE+DELTAY FINAL 
0.318238E+00 0.197901Е-03 0.142380E+01 0.111979E+01 0.338728Е-03 0.142374Е+01 RESULTS 
0.333082E+00 0.280479€-03 0.142414E+01 0.112352Е%01 0.480121Е-053 0.142434Е+01 
0.348089E+00 0.142470E+01 0.112726E+01 0.604827E-03 0.142495Е+01 


0.35328BE-03 


Figure ES. Example Showing the Format of File "outpt3" 
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TRCON = 0.000000E+00 


THZERO = 0.100000E-04 


PO = 0.726000E+02 


TO = 0.820000E+03 


TWALL = 0.520000E+03 


XTEND = 0.265000E+02 


DL = 0.625000E-01 


GAMMA = 0.140000E+01 


OMEGA = 0.76000E+00 


WATE = 0.290000E+02 


CP = 0.240000E+00 


PR = 0.720000E+00 


XINC = 0.000000E+00 


2. Flow Conditions and Nozzle Dimensions 


(Reynolds number value or "x" length 
depending on the option set in the 
Transition Indication Mode) 
(Momentum thickness at nozzle throat) 
(Stagnation pressure) 

(Stagnation temperature) 

(Wall temperature when INDTW = 1) 
(Extension length past nozzle exit plane) 
(Increment in "x" of extension length) 

(Ratio of specific heats) 

(Viscosity- Temperature exponent) 
(Molecular weight of gas) 

(Specific heat of gas) 

(Prandtl number) 


(Not applicable; hardwired value) 


. Input Data from "outpt2" - Data generated from the inviscid contour 


software "char6". 


. Output Code - Constitutes the flow properties along the contour of the 
nozzle calculated by "nbl6" for each set of coordinates that define the nozzle 


contour. 
X (IN) X - coordinate of contour point 
S (IN) Contour length 
Y (IN) y - coordinate of contour point 
M Mach number along contour 
REL (1/FT) REL = (pUg)/u 
RETH Reynolds Number; RETH = (REL)(THETA) 
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PE/PO Static Pressure / Total Pressure 


UE(FT/SEC) Velocity 

MUE (PSEC/FT2) Absolute Viscosity (1) 

RHOE (PSEC2/FT4) Density (p) 

TE (R) Static Temperature 

TW (R) Wall Temperature 

TAD (R) Adiabatic Wall Temperature 

TE/TO Static Temperature / Total Temperature 

DELTA (IN) Boundary Layer Thickness (0) 

THETA (IN) Momentum Thickness (0) 

DELSTAR (IN) Boundary Layer Displacement (8°) 

H Form Factor (H=8*/0 ) 

TW/TO Wall Temperature / Total Temperature 

TW/TAD Wall Temperature / Adiabatic Wall Temperature 

TAD/TO Adiabatic Wall Temperature / Total Temperature 

DPDS (P/FT2) Change in pressure with respect to surface 
distance 

DMDS (1/FT) Change in Mach number with respect to surface 
distance 

DUDS (1/SEC) Change in velocity with respect to surface 

distance 

DTHDS Change in momentum thickness with respect to 
surface distance 

TAUWAL (P/FT2) Wall shear stress (Tw) 

CF Friction Coefficient 

Q (BTU/FT2/SEC) Local rate of heat transfer 


HA (BTU/FT2-SEC-R) Convection heat transfer coefficient 


5. Final Results - Modified contour output as calculated by "nbl6” 


X (INCHES) x - coordinate of modified contour 

DELSTAR (INCHES) Boundary layer displacement thickness 

YCORE+DELSTAR Modified contour accounting for boundary 
layer growth on contoured wall only 

MWALL Mach number at the wall 

DELTAY Modified displacement thickness as 


determined by the equation modeled by 
Demo [Ref. 28: p.24]. 
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YCORE+DELTA Y Modified contour accounting for the 
boundary layer effects contributed by the 
side walls 


C. PROGRAM LISTINGS 
Listings are given as follows: 
1. Program "charó" 


2. Program "nbl6" 
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Program "charó" (cont.) 
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JI ON] 
(TW) LVWYOJ 
Н31411, (2117.)0ч34 
(,NOILO3S 1S31-Q 'LIX3-l,'//'X1) ANWuOJ 
(69 '.) 3LIUM 
N3HL (€ ^O3' NOONI) JI 
JI QN3 
106 01 09 
N3HL (,2, '31^ L9HdAL “ОНЧ” ¿Ye "397 198Н45Х1) 41 
JI QN3 
(06 01 09 
N3HL (+Z: '3N^ L9HdAL *ONV” ,1, ^3N^ LOHdAL) JI 
ЗІ она 
(IV) LVWYOJ 
LOHJAL (8b*.)0V38 
(,1408Н1-2 “11Х3-1,79// “Х1) 1УН804 
(Lb) ILIUM 
N3HL (Z '03* NOONI) 41 
41 ON3 
906 OL 09 
МЭНІ (,2, 7317 20N4AL “Онч! ,Ч,. "39" 2084Х1) 31 
JI анз 
906 01 09 
N3HL (,.6, 
"змо ZONdAL “UNV" Z; ‘IN ZONdAL "ОМ" ,,1, "3М" ZONGAL) ЗІ 
41 ON3 
(1М)1МН804 
ZONAL (9»*.+.)0v3Y 
t,1NILN3NOdX3-€ *YYINDIV-2 *LIOHS-1,*//*X1)1VWWOJ 
(Gb '.)3LIUM 
N3HL (1 'O3* NOONI) JI 
(TI) LYWYOg 
NODNI (bb *+.)0V3Y 
(,АТМО 518104 ТТҮМ”Ч 74315-14-4315"Ч , 
21014100 40 34Х1-У,“ХЄ7/",6508НОНЧ 153174 ’LIXI'V :оч3ч, 
438 01 Н19831-6,%ХС7/”,1Ч08Н178 711Х327У :OVaU 3A OL LHOIIH-Z, 
'XC'/',1NILN3NO4X3^2 ”ЧҮ10938"Ч "'їчон5'ч :31220Н 40 343Х1-1, 
"“Хс”у ле 21) GZONWHD 38 OL INTVA INdNI,'//'X1) LVWNOJ 
(Cb ’o) ALIUM 


NIKL (Z: °O3° ISUNI) JI 
JI ON3 


$06 OL 09 


N3HL (,7, 737: 44581 "ON? ¿Veo "390" JUSNI) JI 


41 


506 01 09 


ON3 


N3AL (,Z2. ^3N^ dSUNI “Онч! ,1,. ^3N* dSuNI) JI 
(1%) 1ҰНЧОЧЗ 

JSUNI (Zb'.)Q0N3U 

'S3A-l 'i123Hu402 S3Q1VA LRdNI 34v, '//'XxUl) LNWOJ 
(6€'») 3LIUM 


AUWYSSIDIN JI LNINI OL SIONVHO 


31 0N3 
(ТТМ 91477 “Х1) 1мнн04 
244100”, - 1МИНО4 104100, (8901's) 3LIUM 


3513 
(21М791М7/ “ХІЗ1УИНО4 
ILdLNO’, = LYWHOd LNdLNO, (L90T’ so) SLIUM 
N3HL (,t, 303" INOJAL) JI 


JI QN3 

(LIN 'O2V '/ '* XU) LNWH0J 

ZHLOT’, = GOvV3u 38 OL HIONTI, (9901'.)31I'M 
3S 13 

(ртм огч ИХТ) миоз 

15191”, QV3u 38 OL HI9N31, (S901'.)31ItM 
N3HL (,1, 7937 Н1144Л1) ЗІ 


JI QN3 

(9v '02V '/'XU) 1v WmioJj 

Qv3d 38 OL LHOI3H, ($90l1'») 3LIUM 
3s13 

(pw‘ozw’/‘X1T) Lywuodg 

ОЧЗЧ 3d OL LHDISH, (€900's) 3LIYM 
нзні (,1, 3037" 19н4А1) ЗІ 


JI анз 

(114 6% *//*X1) LVWYOJ 

6З4Д1”, - 317208, (2901".)311чм 
35713 

(LN '62'//'Xl) LVWHiOJ 

21411", - 31220N, (I1901'.)31IUM 
N3HL (,Z, 703" 208411) 41 3513 

(SV ’6W //'’XT) LVWNOI 

13dAL’, = 37220N+ (0901’.)3LIUM 
NIHL (,1, “037 ZONdAL) JI 


SLINSIY LNAINI AWIISIO 


31 анз 
bo6 OL 09 
N3HL (,2, "317: 100441 “0897 ¿Mo "39" 100411) 41 
41 анз 
v06 OL 09 
NIHL (.2, "зм? LNOJAL “ONV* ¿ls “IN” LNOJAL) JI 
(IW) LvWuod 86 
1nodAL (gt'.)Qv3u 
(,А1НО 518104 11ЧМ-2 41 
,/“4315-314-4315-1 “104100 31М210М1,7// “Х1) АМН80О4 (6 
(LE'«)3LIUM b06 


зі Ома 
£06 OL 09 
МЭНІ (,2, 731" HITIAL ‘ONY’ гм, 3329: ні14А1) ЗІ 
31 анз 
соб OL 09 
мані 6,2, "3N° HLIdAL ONY’ l, “зм' Н11454) 41 
(19) LVW4OJ 9€ 
HLIdAL (96 7.) ОЧ3Ч 
(,/5088НОНЧ 1531 40 18Ч15-2,7/",“Н1їОН31 11Х3-1 “0438 4 
‚38 OL HLONIT TVLNOZIUOH DILSIVILIVIYHD ALVIIONI +*//*X1)LVWYOJ SE 
(SE ’.)3LIUM £06 
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Program "char6" (cont.) 


{S’013) IWW0J 

XIHLA (YONUIFLVLSOL *19*.)0V3Y 

(.453НЭМ1) INOUHL JO LHOTIH JTYH Ілан, ‘и "х1)1чичоз 
(29 ' ,) 3LIUM 

0230393 


JI ON3 
026 01 09 
NIHL 10 7197 ЧОЧЧЗ) jl 
(S'0Ol3) LNWu0J 
X3HLA (ЧО08НЧ3-141501719",.)0ч34 
(,(S3H2NI) LIX3 JO LHOIIH 4J1WH LAOdNI ,'//'xl) LNyniOÀ 
(665 5») 311348 
0-40443 
N3HL (1 303" НІХОНІ) 41 


31 0N3 
JI ON3 
616 0109 
N3HL (0 719% HONYI) 41 
(с"013) 1мичоз 
LXX3X (YOYUI=LVLSOI 95 *+)0V38 
(,153НОМ1) 508НОНЫН 1531 JO DNINNIOIB LNAINI.+*//*x1) LVWYO4 
(15 7.) 31148 
0zuOuu3 


31 ON3 
015 01 09 
N3HL (0 7197 YOouy3) Jdi 
(57013) 14084404 
1ХХЗХ (408834141501796 ,»)ОЧ38 
(+ ISIHINI) HLONIT LIX ININI.'//'XT) LUWNOS 
(6557.)31148 
O=UOUY3 
N3HL (1 *03* ILXONI) ЗІ 


JI ON3 
506 0109 
N3HL (0 "19" YOYY3) JI 
(S°O13) Lywu0g 
1845 (YOUNI=LVLSOI ‘pS’. OWS" 
(,438Wü0N HOWN NOILOIS LS3L LNANI.'//‘X1) Lypwuod 
(£6',)3LIUM 
КЕТУГЕ 


jl ON3 
005 OL O9 
МЗНЇ (0 7192" 40443) ЗІ 
(65013) 1УМЧО4 
ҰКХҰО (ЧОоЧНЗ2: 191501 ‘2с ',) 0ҰЗН 
(,SILN3H 21310345 3O OILYW і1пані, "/// "х1) 1мумоз 
(15*+)3LIYM 
O=-Y0443 


(46099. 5Ч1414М1 el 

.3S3HL 803 518104 TMII30 BZUTILN 2513718 +... .*X9*//”X1) LVWYOJ 
(22*+)ILIVH 

(¿ITZZON 3HL 403 ¿2 

¿«SNOISNIWNIO ANLIHOID ¿*X91*/* .-DIILDIAS ONY SIILSIVILIVIVHI +1 
38013 44312346 6104МІ 30 135 ЇХЗМ ЗНІ,'ХФ"////''х1)1мичО3 
(127-)311848 


SINANI AYLIIHOID DIJIDIAS ONVY 52115143104Ч8ҮНЭ 8013 31220Н 


25 


15 
005 


22 


12 


2=LNOONI 
3513 
T=LNOONI 
нн *03* LnOdAL) 41 


JI ONJ 

C* 3LXONI 
3513 

1531ХОМ1 
“O3° HLTdAL) JI 


31 ON3 

C*HLAONI 
3513 

I=HLAONI 
703" 198Н4Х1) 41 


31 4н3 
£-uSONI 

3513 
Z=WSONI 

“037 208411) 41 3513 
145081 

Col. 703° ZONAL) 41 


u3931NI OL NOISU3ANOOD 


31 083 
0601 01 09 
NIHL (1 *03* 100М1) ЗІ 
31 ОМЗ 
0201 01 09 
N3HL (Z '3N' 100М1 'ONW' T '3N* LDONI) I 
(11) LUWYOJ 
LNONT (116 *.)0V238 
(,ОН-2 *SIA-1 “С539МҮНЭ u3HL10 ANV, '// 'XTl) LVWuO. 
(016'.) 3LIUM 


ЗІ ОМУ 

(LIW‘STV’/ ’XT) Lywuog 

214100”, 104100, (8:01".13114М 
2513 

(21ч 914 79/ “Х1) 14н4034 

114100”, 104100, (440175331188 
мані (,(, "032" 100444). 41 


31 UNI 

(41470247 / ХІ) Іо 

28191”, - 0438 38 HLONIT, (9(017.)3118М 
3513 


(01Ч70247/ “Х1) 1чуг404 
1Н191%, е 0434 386 HLON3T, (SLOl'.)3L1uM 
МЭНІ (,1, "03" HITIAL) 31 


31 0N3 

(9v’ozw'/‘x1) Lwwuod 

ON3u 38 1НЭ913Н, (bLOl's) 3LIUM 
3513 

(ҮҮ70247/ “Х1) 4404 

OVIY 38 1Н913Н, (Є(017.)311848Я 
N3HL (.1, 703: 19Н411) 41 


41 анз 

(LIV*6N *//*x1) LVWYOJ 

(34117, - 31270М, (2:017.131148 
3513 

UN *6N *// RU) LVWIOS 

234117, = BIZZON. (1:017-)3118Ы 
NaHL (,Z, "03" 208411) 31 3513 
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f, l '*ON OuN2, 'X02'//*' ,М1ЧО 10481, “Х02 )i1MWUOJ3 Ib N3HL (€ 7037 DHINI) JI 
JI ON3 
10hHT130 '12u 'X3HLA ' LXX3X ' LHS 'VWWNODT 
“аЯ20М1 “4У10М1 'TVOQNI '1n00NI '31XONI 'HLAONI "HSONI (19: 791) З118М JI ONG 
(////9Ч0Т7”Х02) 1М8Ч03 906 StS OL 09 

(01 *T=1 (1) 1H3) (909 *9T) 3LIUM N3HL (O0 '19* чочча) JI 

(6*0TI3) LVWWOJ 

185 (408НШ3-14150176(7.)0ч34 

(GAWS ЧЯ04 30ЛУЛ MIN LNANI,*//*XT)L1VWyJO4 

(8 *+)ILIVM 

0zuOuu3 


(/////*' ,317220N 3HL 30 3402,€ 

JSIJdOuLN3SI 3HL,'XO0Z2'//',I1 Luwd,'X*t'//',S931220N OINOSWM34fS ,Z 

‚IVNOISNIWIO-OML 30 НЭ1530 IHL dOj4 WWvu5Oud ua3lündHOO, 'xG'//'I 
¿ONVIAYWH “MYO ILIHM ‘AYOLYHOEGYT JONVGYO TWAYN "S'D, 'XST1'XTI) IVviO3 Ob 


Program "char6" (cont.) 


(0$ '9T) 3LIuM 


JI ON3 
JI ON3 
777 108 O1 09 
мані (T 303" LX3NNI) JI 
JI ON3 
208 01 09 
NIHL (Z “IN” LXINNI ONN’ | °3N° LXINNI) ЗІ 
(11) LVNyYOJ 
LXINNI (26*+.)0VINY 
(,0N-Z 'S3A-[ 'LS39NNHO VIHLO ANV+*//*X1) LVWIOJ 
(16 ".)З1Гїчм 


(1/*5"013*6W*x1)5"/'5"013*6W*XxI *//*xT)LVWdOJ4 
10130”, з 1nWwI3O0,'l10u', = LON, XIHLA*, = XIHLA. 
“ҮХХЗХ”, < 1ХХЗХ, 166%, з 1Ғ5,”ҰНИҰ9”, = чуучо, (96!) ILIUM 


31 ON3 
JI ON3 
696 0109 
NIHL (0 719" ЧОНЫЗ) ЗІ 
(5* 013) LVWYOJ 
1nWH130 (udOMd3-1VISOI'68'.)GX3W 
(,10W130 4034 3n1VA мам 1п4мІ,"// "хт) 1мичоз 
(88",)З11чм 
0O=YOWY3 
NIHAL (9 *03* OHONI) JI 


JI ON3 
31 анз 
096 01 09 
N3HL (0 "1907 ЧОЧЧЗ) 31 
(5*013) 1WOJ 
LOY (ЧОННЗ-14Ч15017(8”,)0ч34 
(¿LN YO ЗПЛУлЛ MIN IQdNI,'//'X1) LMWuOJ 
(98'.)3LIuM 
0zdOWu3 
N3HL (g 'O3* онон!) 41 


JI ОМЗ 
JI QN3 
566 01 09 
NAHL (0 7197 ЧОНЧЗ) ЗІ 
(5*013) LVWy4OJ 
XIHLA (YOYYI=LVLSOI *p8*.)0VIY 
(¿XIHLA моз INTIVA MIN INANI¿*//*X1) 1VWyJOJ 
(368 5.) З118М 
0-uOuu3 
N3HL (b 'O3* 5HONI) ЗІ 


JI ON3 
JI QN3 
066 01 09 
мэні (0 "19" ночиз) ЗІ 
(5 012) 1мичоз 
LXXIX (40843-141501728”,)ОЧ34 
(LXXIX YOJ INTIVA MIN 1п4мІ,"/ / "х1) 1мичоз 
(18) ILIVM 
о=чочиз 


N3HL (£ '03' 9HONI) I 


JI ON3 

31 ang 
095 OL 09 

N3HL (0 319" yodda) ЗІ 
ай (5"013) 148404 
VWWVWO (МОНН3-1Ч15017:(".»)0ч38 
(,VWWV9 HOJ INIVA MAN LNANI,*//*X1)1VWyOJ3 
(9: 4) 3114м 
"ЕКСЕ 

N3HL (1 03" 9НЭМ1) I 

(11) МУЧО Si 
ЭНЭМІ (5:".)Омач 
(,108130-9 “1084-4 'XIHLA=-b "їххах-є , { 


(“145-2 "чинчОо-| :азончно 38 OL INIVA LNANI,*//*X1)1IVWYVO3 be 


(be “+)ILIYM 108 
NIHL (Z *03* SNVNI) JI 


41 Она 
008 OL 09 
N3HL (Z '3N' SNVNI “ONV” I °3N° SNWNI) JI 
(TI)LVWYOJ £L 
SNYNI (Е('.)ОЧЗЧ 
(,0М-2 *SIA-1 *¿LDIVYOD SANTVA LNANI IV +*//*X1)LVWUOS Zi 


(Z¿*+)ILIUM 008 
AUWSSIIIN 41 SINTVA 10481 ЗЭМҮНЭ 


(1(/“57013“64?Х1)67/7“6:01376У“Х179// “Х1) 1УУчО4 11 
108130”, - 104130,'124', - 1028, 'ХЗНІА", « ХЗНІА,І 
'LXX3X', 5» LXX3X,'LWS*', - LWS,'VWWND', = VHWVDO, ((( 7.) 3118М 


811053Ч 10481 AW'dS!ü 


Зі ма 
565 OL 09 

N3HL (0 “19” YOYYI) 31 
(S" 013) 1VYAYOJ 89 

TINATIO (YONBUI=LVISOI “89 *.)0VI8 

(,(514 009-4"0 "514 051-031) .І 
«LNIHIVINI X ноз WOLlOV4 JAILVIITALLTAN LüdNI,'//'XI) IMWuOA L9 
((9'.)ЗіІЧм SES 

Озчочиз 


31 083 
066 0109 

N3HL (0 7197 ЧОНЧЗ) 31 
(6013) 1ММЧ03 99 

19H (uOuu3-1Nl1SOI'99'.)0V3W 
(,(S3H2NI) LVOYHL 1Ч ЗНПЇҮЛЧПЭ JO SüIOVu LüdNI,'//'XI) 1VWuOJ v9 
(99 “»-) З11ЧМ OES 

0-404Ч3 


31 ON3 
ЗІ ON3 
925 OL 09 
NaHL (0 "45" нчоннЗ) ЗІ 
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(+ON=CN 
or*t=f 066 00 
ON-ONIN=OL 
(QN3N "1HS 'IX 'ON'MX)OQ3XLI3S 1142 
0” 0” (ONIN)MIHL 
LWS = {ON3N) MWS 
“14 (0838) МА 
(QN3N) 1X» (ONGN) MX 
JONILNOD 
(ч10172) ч13а (чЧ10171) ч13 
(VLOI '2)3HL= (WLOT *1)3HL 
(ч10172) 55 (ч10171) 85 
(9101 72) Із (9101 '1).1 
(9101 *'Z) Vs (NOI 'U V 
(VLOL ’ZIA= (VLOL’T)A 
(М10172)Хс(М10171)х 

OX*I=WLOT 006 00 0S8 

(9"S13b 'SIZ*+.H1*X1)LUNYVOI ZEF 
OHL'(N)MHS 'AA'XX'3'N (26b '90) ILISM 
£^LS« IN) MAHLEQHI 
LIX34« (N)MA2AA 
LIXIA +» (N)MX=XX 


OLb OL OD [T° LON) AI 

("T-Zeo(N) INS) LYOS/*T={T‘Z)¥ 

(3) 7Ха(172)Х 

(4) 15454 (172) н5 

046 01 09 (М3717М45Ұ7) 41 

N= IN 

N3N‘T=N 006 00 
(/",Зпому нчачісз,'хс','ом нами, "хі", (53НОНІ)А,'хс"2 
« (SIHINI)X,*X9 * 5, XB ON, XS '// 4 (SINTVA TIVM OILVIOÍBILNIT +) el 
: NOILYINdINVH SIILSIVILIVUVHD.,*X01 *///*Xx1) Lvws03 111 
(11191) 3L 19M 

0*0=(1'Z)v13 
0'0=(1*2)1 
0*0-(1'2)3HL 

0"0-1172)3 
1-ONIN=NIN 

гзчаачи 
ОН 61 094 NOISSVY OLI 
(108130 “мч13” 


LIXIX 'VHWVO “ONIN “NISV1*8 “LIV LX *LIXGA НІА “LHS “чОЭ2Х “185 “1Х) ОН5 TIVI OO! 


011 0109 
(10нм130 “мУ13 1 


(я "(мумані "(2-4 "2)ч13 "(11 LIX3X “VHWNVWO ONAN 'NLSVT *8 "198 “LX 'L11X3A “HLA *LHS “YODX *1HS “IX)OAX TIVI 86 
-1*2)W13* 14*2)V13" (2-42) 3H1* (1-M'*2)3H1*(%'2)3H1'HVW13) YILNIM TIVI 011 01 09 
(451) мн57 (22-32) Ч11 | (108130 “мм13 71 
3 (т-х 2) 443 “5(Я“2) 137 (2-4Я“23)н5“5(1-Я7"2)н6 5 (м72)15 “мү13) ЧЗ31М1Мм 1192 11ХЗХ “ҮнННЧ9” ОНЗМ “м15У178 “128 “1Х711Х34“Н1Х 7185 “ЧО2Х “115 “ЛХ) 9384 1142 S6 
(WON) MAS (7-H гт B6 OL 09 
уміЗ "(1-м "2)мІ13 9 (572) ч1 39 (2-я“2) Х51(1-я 2) АЛ (2) Х.мм13) ЧЗ1М1М 1142 091 (,М3Ал19,2Х17/7,104м1 ОЗ 314 1440444М ;1 
IN=ON «LON Si NOILAYIVLSIO "ON HIVH INITIBILNID DIBND+.+. XL) LWWHOS 0021 


Program "char6" (cont.) 


NH OL 094 NOISSV 
(092 ‘ObL) ‘AW OL оо 

(ХМ) мх (г-и гг) Т 

vV13'(U-43'2) N13' (4 2) N13' (2-43 '2) X' CU-43 2) Xx " M' 202 X'MNI3) u3LNIM TIVI 
058 01 09 

INNILNOD 

024 01 09 (мХ13'19" (Х'2)Ч13)4І 

(9“51397512 хі) ічнчоз 

oHL'(S'Z)HS'AA'XX'W4'N (169 *91) 3LIUM 

€" ¿S- (47) 3H L=0HL 

LIXIA + (N'2)A=AA 

11Х3ХА.(472) Х-хХ 

059 01 09 (27037100081) JI 

2“Ч444ЧХ-ХЯа414 

(0472) ч13 “(57235 11-32) 6 (1-3Я72)31 

(A ZIJA  CL-23'2) X" t0 '20 X" Ut-M'2)NI3' (n 'z)3H1L'05'Qvws 'üaon'v2) 3a1V 11v2 
(11-35 72)Х-14 52) Х) «142411-472) А41472) 3 

(LWO-1 WD) / (01-41) XeLVO- (1-1 2) XL WD* (1-X SZ) A= CTH TAD 2 (1 '2)x 
(УМ о VYL-" 1) / [VVWVW+VYL) -1у0 

(89. ву ўе" () / (вуў-891]: 172 

“2/ (1572) 54 (1-3 52) 6) 6995 

(я52)М/ (5152... IR IV) LUOS=- IN ’ZIMS 

2/4032) ме (1-3 71) У) МУҰ 

"2/460 '2) 1* CUI 0 D) 2VVL 

%2/((47%09)945(1-Х72) У) <8УУ 

“4/ (1372) 19 11-472)1) 54ч1 

(1572) ч “Ч1110172543 “Ч3н113Н1 “0мМ7(1-372)Ч 11-34 1) уу ноэч TIYO 
(1872)3Н1) 50Э/ ((Я72) 3нН1) М15-1472)1 

*2/(43HL>73HJ) s (4*2) 3HL 

((t-23'20N/* UNVXN- Ct (E73 '2) Vs0M) / U NNIN eMe CL-3 2) 3HL2 83HA 
((t-43'0)V/^ONVINe CELO D ven) /^ t NN1N enM- (CC D 3HLST3HL 
Язох 

н'чадчівя 00: 00 

l* vdd Vx sv da voi 

06b OL O9 (N'LO'NLISVT)JI 

(vdd V3 '2) A, MV13- (NVdd X ' 2) N13 

0°O= (Wdd¥X г) зні 

0'Оз (маамн'2)1 

VNVIs (dd NX '2) V 

IWNS= (vaam ’ZIns 

YVNYL. (IX INITXI = (VAN 2) A 

{N) 1X= (wddWu'Z)X 

058 01 09 


123 


924 
004 


059 
16% 


06» 
04% 


ott 


(0021 ’.)3LIUM 
(9'2013', 9 'мзлі9 Іпамі чоз зічІчаочаач 10м 5І1,'ХЕСІ 
*/* ¡NOILNBIVISIO “ON номи INITUILNID JINI »..*X61 /// *X1)LVWHJOS 06 
g (06*91) Зї18м 
66 01 09 (07171178 0а8У7772-71978)31 
0" ¿- (IHLA/LIY/MD) LYOS + (071-186) /072/1Х:8 68 
Ч50М1 (867587001) 01 09 
ЗПМІ1МО> 08 
MVY13=(1)MA 
00000°t=(1) MWS 
00-11) мХ 
009*1=1 08 оо 
10“1-185-.115)1805-11Х3А-11Х3Х-1Х 
1ХХЗХ-11Х3Х 
06 01 09 
10“1-1445-185) 1805,11Х3441Х:11Х3Х 
1ХХ3Х-1Х 
69 0109 (1703! ЗІХОМІ) 31 
MVL3 +XIHLA=HLA 
XIHLA=L1X3A 
09 01 09 
XIHLA=HLA 
МҮ13/Х3Н1Х-11Х3А 
05 01 09 (1703 “НїХОМ1)41 
(°[-Z. eLWS) LYOS/ * T=VNVL 
WYUILeLHS=MV LI 
0Dr. (Ze. LAS :«89*0"Z)/VO)=W83L 
(39) LYOS=NM 
"2/ч9з19 
89/%2-99 
89/ҰҮ9-39 
Э9/Ч9-09 
49.“ 2«-29 
“Т-ЧнННчО0.89 
"Течуумогчо 
62аЧ11101 
t00* 022543 


(S'213',» aünwi130, 'X»£6 

!J!6* 213,5. 124, 'XtC"/'S" 203, x3BliA,'xcc'/'S"2z03',» іххах,9 

“кб '/'S* 213',-. iNS,'XLC'/'S^ 213^ ,»* WANWO, ‘XS ",2 “ОМ OUVI,'X0ZS 

//2X1',* Ou2QONI, 'Xb£ '/'211',» aMLONI , *Xb£8 

*'/2VI ';s TNOONI , 'Xb£ '/'201/^,s LRnOQNI , 'XbC '/'2T1I',» SAXONI , "XP£L 
'/'2V1',5 HLAONI , 'Xp£ '/'2U1',7 Ч5041, "Хб 





Program "charó" (cont.) 


£8 OL 09 
NN OL €8 NOISSY 
6-38615670-492 
2-321529"08562 
11410239 2112) 
152:170462 
16294704422 
Рс6%-Тг12 
28 0109 (%110971731722Х) 31 
MWL3/ (I) 1X2D9K 18 
(68718) “АМ OL 09 
X130+ (1-1) 1X= (1)1X 08 
10H130.x1302x130 
LIX3X»^09b/4(^ U* ULIX3X/(U- I) 1X) 4 000) 5019x130 
08 01 090 (02°371°1) 31 
1441 
0097241 88 00 
1000: 04х134 
0"6-11Х3Х 
LIXIA+8X= (1) YODX 
8NS= (1) IWS 
8X- (1) IX 09 
3NNILNOD 05 
09 OL O9 ((000*1'39'8WS) JI 
(2X) 9414s «gus 
МЧ13/Яхаох 
10000*0+8X=8Xx 
009 "Тасі 06 00 
0*0=8X 
AN OL 18 МЭ155ч4 
2-3%)98170-:92 
(68600452 
60tbb*0-*p2 
£8008*02€2 
269620 0422 
071412 
(2%Х.(22%Х.(62%Х» (Р2%Х65(62%Х.92))))) г (Х) 93135 
(009) чаэх”? (009) 145 “(009) 1Х NOISNIWIO 
(2-0'H-v) NOISIO3ud 318000 LIOITdWI 
(INTO MY LI ' LIX3X'l 
VWHV9 'ON3N'NISV']'8'L2u 'LX' L1X3A'HLA ^ LWS ^uoOx "Tus "1x) Qus INILNOYENS 
ON3 
Nunil3u 
(9 °613 хе '9 "613 °с̧І'х9т) 1мичаз OOT 
(акак ‘тг * (A) "INS (4) заЭх “) (001791) 3118м 
LWS = (ONIN) IWS 
LIXIA + (ONIN) TX = (ONIN) YOOX 
LIX3X* (ON3N) 1X 
WLO=ONIN 66 
JANILNOD 06 
66 OL 09 (LIXIX *'39 *(()1Xx)31 
LWS# (Cf) INS 
11Х3Х,(Г) 7Х41Г) ЧОЭХ 
XT130+ (1-0) 1X= (0) 1X 
10н134.Х134-х130 
LIXIX + "099/41" T+ (LIXIX/ (1-0) 1X) +. *001)901=Xx"130 
с=с 
009 ’IN=eL 06 00 
IHNISYI=IN 
LIX3A« (IM 1X8 (1) 402X 
Lx= (1) 1x 
LWS= (I) IWS 
LIeNLISV1 $8 
INNILNOD 88 
$8 OL O9 (L1X'39* (I) 1X) JI 
( (1)1x) 9413s * (I) 18S 
LIXIA + (1) 1X= (1) Y0DX 
x130* (l- I) 1X3 (I) 1X 08 
10130, X130«x130 
LIX3X.'095/ (C I* (LIX3X/ (1- I1) 1X) .* 001) 2015X130 
08 OL O9 (0Z2'31'I1)JI 
I*4I 
009'zsI 868 00 


1000  бех130 
11Х3Х,9Ха(1) ЧОЭХ 
8WS» (1) TAS 
8х- (1)1Х 09 
INNILNOD 0$ 
09 04 09 {1000°1’39’aWSs)4I 
(9xX)9JINS=8HS 
10000*0+8X=8X 
009 1451 06 00 
0*0-8X 
LIX3A/LIX3X*LIX3X 
11Х3Х/1Х41Х 
(/' ,"oN HOW, Х617,Х.4х(7.М.4х0279772 
? NOILWNOZT DIGND ONISN NAOILNGIYLSIO YIBGHWAN HIVH INITVILNIO +1 
'x0t'///'c^t€t3z'Xz2'/',8, 'X0l' ,.19u, 'XLC'XTl) 12NWdOÀ OZ 
8g'IOu (OZ '91)3LIUM 
(/'s*tti3e 'xi'/'z 
~ GLX. KO! LIXIXs XL ALIXIA, XOT HLA. XS SON HOWH X6' //'l 
«(SIHINI NI SHLON31) sual3WHvuvd 31220N,'x61'////'XQl) 1vwuoOJ4 O1 
LX'LIX3X'LIX3A'HLA'LWS (OU'90) 3LIUM 
(LX/X-0* 1) «e (LX/X-0' D e (LX/X»8-0* 0 « (0 L-LWS) -1HS» (X) 24193S 
(009) u02X ' (009) 1S * (009) 1X. NOISN3HIQ 
(2-0'H-V) КОІ5ІОЗЧ4 318000 LIDIIdWI 
(108134 “му13 “11Х3Х 71 
VWWV5'ON3N 'NuisV1'8 '12u 'LX' LIX3A 'HLA ' LWs “ноЭх “ме “7х) 938. 3511008805 


ONY 

4015 
(,Zadano 3HYN41IJ NI,'/',Q03uOLS 51 379149 нҮнН008Н4 ,'t 
650026 ІЛ IHL NI O3SN 38 01 3114 51Н1 WOU3 V1V0,'///'Z 
134380 ЗНЧНЗ114 КІ 034015 «1 


г4001НО3 31224МН чаз VIVO DIdOuLNS3SI 10d100,'///'Xxl)4vWdO4À 0011 


(001174) 311%М 
(11) 35012 
(91) 35012 


(1 онїмїч 

LU 3114 ON3 

(0236) Lv WdOÀ 

Corser (Cy IWS ‘LNOD IC) MWS’ (COMA ICD IX) CTO 0U' LO)3ZIuM 
(£°OZ3b SI) LYWY 

LUXAX LX LIXBA ‘HLA ‘OC COUT LUD aALIuM 

(8WOT) LVHYUO З 

(01 'їзг'«суїна) 1291“11)3118мМ 


0°0=1N09 
(9°S139 ‘SI “XI)LVWUOJ 96) 
(00 *121*(1)8* (1)M3H1 * (1) MHS * (1) 15 * (1)MA * (1) 1X*1) (b6b 91) ILIVM 
ү, (NI) Auün2 40 Qvu,'XZ2',319NM TIVM,'XC',"QN HoOMW TIMM.Z 
'X£','ON HOWVAW 12,'xXS ' , (SSHONI)A,'X9' , (SSHONILI X, 'Xb ' ,. I, 'XG*'//* ,t 
uUnOlNOO 3402 OIdOuLN3SI “1 1844 04 104100 1WVNIJ,'XOC '///'X0) LVWHOJ3 66» 

(£6p '9t) 3L TUAM 

0*0=(00)Y 

0 '0= (т) ч 

INNILNOD 

€" £S»IHLO/SOz (PC) Y 

9-3" 1=3HLO0 (0*0*03"*3HL0) 41 

(1-9Г) МЗН1ї-1(1»0Г) м3НА«3Н10 

(264 (11-0Г) мА-1(144Г) МА) 5206 .4((1-ЭГ)МХ-1149Г)Мх)) 1805:-50 

lle ’z=ber [66 оа 

l-or= ttc 

4. 1Х3Х.1Х:1Х 

11Х3Х-11Х3Х-11Х3Х 

3üNILNOO 

C°LSe (CN) M3HL-( с) мэні 

LIX3As(CN)MA-( CIMA 

11Х3Х,(ГмМ)мх-1 С) мх 

1ІХЗА, (см) 1х={ Г)ЛХ 

(сМумибе ( г) мис 

ICN)mSst с) 195 
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Program "char6" (cont.) 


(X 'auaL ^ux'cx'Ix'L's*u'ux3) YILNIM INILNOYANS 
ояз 
нз01394 0611 
(,1мПп0Х,“Х0О2) 1МНЧ04 9911 
(991191) 3LIVM 
0511 01 09 (ч11101!117531аг) 41 
(53102) LONY=2Y SHUI 
Занїїноэ 0111 
01 
111 “с91175»11(2543-((1н003) 10нМ/ ((1-1800) 1ОМНҮ- (18003) 10ММ) ) 544) 31 0011 
INNON=S31NC 6601 


“2/ (3-5) 4((1МП0Х) 10ММ/7 1) ММ1Х-(((1М003) LONV+oM) / * T)NVLY o Me (LNNON) IVY 0601 
((Z-LN0O0X) dY- (1-LNNON) IV) / (1 

(Z-LNNON) LONV- (T-INNOX)LONV) » (1-LNNON) IY- (T-LNNOX) LONV= (LNNON) LONY 0801 

9111071 %6-14н003 0111 00 0401 


72/ (4-5) * CULO LONV/* UT)UNVIV- CC COD) LORV M) /* UUNVIV.Ms (2) aV. 0901 
1“? (І) LONV- (2) LONV OSOT 


“2/ (8-5 * LUCI) LONV/" D)NVIV- CCOCL) LONV M) /* I)NNLVeMs (D) aV OpOl 
*z/ (0*4) (1) 10нч 0601 
(OOL) dv’ (001) LONY NOISN3WIO SOOT 
(2-0“Н-Ч) NOISIDINA ITANOO LIDIIANI 
(ом'чі11101152543'ч'є'м'о'да)ноом 3н11008805 
ON3 
Nunl3u 
(9 стаг'єт'х81)1мичоз 
(tQN3N ' Lex * (0) TAS ' (1) 900X 40 (001 "91)311чм 
LHS= (ONIN) IWS 
11Х3Х, (ONIN) 1X= (ONIN) YOOX 
LIXaX= (ONIN) 1X 
HLC=0ON3IN 
3NNILNOI 
56 01 09 (LIX3X ^39 "(£C)1x) JI 
LWS = (£) TWS 
LIXGA» (£) 1x= (0) YOOX 
X130+ (Too) Tx= (CC) 1x 
1nW130.x130-2x130 
11ХЗХ.“095/ (“19 CLIX3X/(l-DC) 1X) »* 000) 5012x130 
ГЕНГ 
009 “1М-Г 06 00 
INLSVISIN 
LIXIA+ (1) 1X= (1) YODX 
LX= (1) 1X 
IWS= (1) IWS 
LI=NLSW1 
JANILNOD 
$8 OL O9 (1xX'39* (I) TX) JI 
COI)1X) 941385 (1) 18S 
LIXIA + (1) 1X= (1) 909X 
Xx 130+ (1-1) Tx= (1) 1X 
1ünHT130.x130-x130 
LIX3X.*09p/ C Ce CLIX3X/ (1 - I91X) . * 000) 2012 X 130 
08 OL O5 (OZ2'31'I) JI 
1511 
009 ’z=I 88 OO 
1000°0=x730 
11Х34.4Х4(1) чОЭХ 
945 = (1) 195 
ӨХа (1) 1Х 
SüNILNO2 
09 OL O9 (1000'1'39*'8HS) JI 
(8X) 9JINS=8NS 
10000*0+8X=8X 
00914651 05 04 
0°0=9X 
LIXIA/LIXIX=LIXIX 
11Х3Х/1Х-1Х 
0*1-1WS=ZHL 
0* U- (29) 4X3I=IWNL 


{/'s’ON HOWN, ‘XOT’ “Ne ’X€2'//' NOTLVNOS .Z 


«TVILNINOIXI ONISN NOILAGIULSIO u39HON HOVM INITUZINGO +1 
а 0" 6" 136 'Х817/".22,'Х8",1ий0, XB’ , 12d4 , "XCC'XT) 2NWuOJ 8b 
. z2'Xama'iou (89 '91) 3LIuM 
| (ZASE CTIS XS" 
AX ,!X8' 4,LIXSX ,'Xt* 4,LIXSA,'XOU' HIA 4 ХЕ, "ОМ НОМИ, *X6'//'1 
4(6ЗНОМІ МІ SHION31) Su313HvVuVd 31220N,'x61'////'xl) 1INWJOJ tb 
IX'LIXSX'LIXSA'HLA'IWS (Lb 'SU)3LIUM 
£9*22 Sb 
dOLS 
(61”,51МУ003 О4Х ONS NI NOILYYILI,*xXO0Z *X1)LVWYOJ Ob 
1311 (0b*91)3LIVM GE 
06 01 09 
1000"0-29 (00:0:03"29)41 
13.18-19=29 
23=13 
29=19 
S£ OL 09 (05 "19 "YILI)JI 
1934311-3311 
(13-23) / (19-29) = 18 
с̧р OL 09 (10000 "37° (23) 58“) 31 
4721-11-23 
(71-43) /43.-29."2-21 
(29)4x3=33 Of 
072/1-29 
l/'t-'lsl3 
O°2=TL 
100*9=19 
0-3311 
(* t- 1IHS) / LX « XOHQs 1 
(Н1Л/12М4/ (7 1%ҰНЫУО) ) 180656 70-хана 
Ча (ні UC 07 CULX/X-^ D e (LX/X-*1).29)9X3)-"1) o. ZWL=(X)OJINS 
(009) YODx ’ (009) INS ‘ (009) 1X NOISN3NIG 
(2-0“Н-Ч) нО1512384 314000 11211441 
(1048130 “МУ13711Х3Х71 


VWHND 'ON3N 'NLSV'T'8 "Lou "LX 'LIX3A 'HLA ' LWS “YODX “INS “1xX)OJX INILNOWANS 


ON3 
Nunl3u 
(96132761 ”Х9ФІРІМЕНОЗ 001 
(ONIN ‘T= * (A) THS (3) чОЭхХ “1 (001791) 3118М 
(7° s°ON нэмн, ХОТ“ хо ХЕМ, "Ог “ИБ 


',3i1220N l8OHS Y YO3 NOILNBIULSIO YIAWNN HOVW BNITYSLNID, XOL'///E 


“6"Е136“Х171",1Х,2 

'X8', LIX3X ,'XL' ,LIX3A , 'XOU',HLA?'XS' ,^ON HOVW, 'X65'//'! 
¿«(SIHINI NI SHLONIT) SualaWvuvd 31220N,'X61'////'Xl) LVWWHOJ 86 

(NISV1)8OOX'(OQN3N) NOOX 'LIXSA 'HLA'LNS (86 *91)3LIYM 

LWS = (ANIN) TWS 

LIX3A« (ON3N) 1X « (QN3N) 802X 

LIXIX= (0ONIN) TX 

WLC=ANIN 

INNILNOD 

56 01 O9 (LIX3X *39 “(Г)1х)41 

LWS= (C) IWS 

LIXIA» (0) 1X= (1) YODX 

XT130+ (1-0) 1X=(0)71X 

10H130. x1302x130 

LIXIX+*099/(* 1+ (LIXIX/ (1-0) 1) ."001)907=Xx130 

CALC 

009 “1М-Г 06 00 

U+NLSYI=IN 

(:1-2»..48М5) 140541Х-11Х3Х 

11Х3Х,1Х-(1) ЧОЭх 

(1) 1х-1Х 

LWS= (1) IWS 

11-М16Ұ1 

JINNILNOD 

68 OL 09 (1HS*3I9* (1) INS) AI 

(9X)93INS=(1) IWS 

(MVI3/(1)1X)90122X 

ІІХЗА, (1) 71х= (1) ООХ 

98 OL 09 

(9x)931INS= (1) HS 

MYLI/ LI) IX=IX 
LIX3A» {I} IX= (I)UOIx 28 
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Program "charó" (cont.) 


SüNILNOO 
Ий5-2%АДачАДО 
12%2%ЛАзАД 
£23/ (I) A22 
SJANILNAS 
£E9+WHNS=HNS 
INNILNOD 
(1Я)Х-ХХ) «25-62 
02 01 09 (С "03" Я "НО" І 2037 Я) JI 
5144 1-4 02 08 
“1262 
oc 01100801 01: fen d 
sLdN'Is^ Ot OQ 
а= н 
INNILNOD 
(1С) Х-(1у) Х) «22522 
(11) Х-ХХ) «13-13 
01 01769 (т "03% г) I 
514М"1-Г 01 008 
“15-42 
*1=19 
SLAN‘T=I 05 00 
0-АА0О 
“ОзАА 
(1DA*(1)X NOISNIWIO 
(2-O*H-W) NOISIDINA ITINOO LIDITAMI 
(AAQ'AA'XX'SLdN'A'X) LNIOVTI INTLINOYNANS 
QN3 
NUnla3u 
INNILNOI 
(CH NNS s (C) IWS 
(P) MXe (C) 1X 
онам'омег 09 00 
INMNILNOD 
(030 “(1)4045 “(1) мх”с5 (9-4) 185 7 (9-3) 13) 18519417 1140 
or OL 09 
(040* (1) HS” (1)MX*S* (2-4) UNS * (2-N) 1X) LNIOVWT TIVI 
ас 9:4 (09 
02 01 09 (ONIN*IT*N) II 
S OL 09 
Y «s 
Ol OL O9 ((3)1X^31' (I) 4X) AI 
ONIN ’ON=I 0% 00 
£-ONIN=ONYM 
LO=A 
INNILNOD 
11 OL O9 ( (ON)MX'39* (C) TX) JI 
CEEC 
009 ’I=C OL 00 
(009) WWS (009) 75” (009) 71x* (009)MX NOISNIWIO 
{2-O’H- ч! Ма15123844 318100 11317141 
(ONIN “INS “TX “ON *MX)OIXLIS 3М11008805 
ON3 
NUNLAY 
(LAZA MATA +? (XZX 0 (XZX) LOS UNS o LIZWSI JOA + ((15) 404) 627 ІМІЗг2УМіЗ 
(285/"1915/"1).4"-мм5 
TA-ZA= TAZA 
[X-2ZX=1XZX 
89.4 (г. РІЗ» 909» *£)/N9) РК (13) 304 
(2-07Ң-М) МОІ5ІОЗЧА 318000 11217401 
(2М137206 %1Ы641А%2А7(Х7%2Х71М137ЧМ13Н1789"8УИ5 "8094У99)31М INI LNOUGNS 
амз 
NyUN LIA 
(1Х-0Х) /54(1Х-НХЗ31)»,»(СХ-1Х1/3.«(СХ-НХ3) “4831 
51 01 09 
С(ЯХ'СХх'ІІ 
х”1757н) 3Х18л10,.(Гх “Нх3) 34104 (Ссх “1хХ7555н) ЗАЛІ), (Іх нх) 33І8гч-4ч31 
01 01 O9 (2 303" Я) ЗІ 
(Ax * 1x1 
“tax 'CX'L'S) 3AVAIQ' (CX' IX'S 'u) JAVAIQ) 3ANAIQs (4X 'CX'IX'L'S*u) J3AGAIG 
(CX'IX) A3IQ/ (S ^u) 33IQs (CX' IX 'S u) JAVAIQ 
CX-IX& (CX ' Ix) 34IQ 
(2-0 'H-W) МО1512384 318048 11217441 


SI 
01 
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Program "nbl6" 


Z=OUdONI 
(442041 
0" 0=3NIX 


а5пм4 


eLNIdd 
«ІМІЧА 
eo LNIUd 
«іміз 


(¿312204 343 JO 10014002 зЧу 59,7/ 75 ,,119М 69 әлетра ,р 

‚apis 342 103 192487 AÁzepunoq 243 YI0Q IUNOIIV OJUI] ,”/”,әхез , 

‚,оз ретітрош чәәд sey ue1boid эці “зэХер 416ёрипо4 943 10 «3э9)19,2 

г, ЭЧЭ, /  зчлооо OIUJ өх Ооу І 3194 шо1) 18013002 Э|ү4013496ү,| 
*, эцэ єзуртрош шезбола 6үч1,9хХ(7// “, шоүзэвролх3и1,7/ “Х1) 14804 61 

(521 7.) 3114 

(,51233443 500251Л 403 0313100н 80001802 38023 2140810351, 

"Х5177/5.11 1МЧ4,“Х9С7//7,31220Н Э1НО583405 ТҮМО15М3Н14-0М1 Ч ,2 

‚30 NOISIO 3Н1 403 НҮЧЭ084 YILNAINOD ¿XS *//* ¿ONVIAUVH ‘NWO ‚I 
«I3LIHM “ANOLVVO8WT зоммачо тмлмн 'S'n, 'XGL'X1'//////'XV) INWWOA OZl 

(02175) 31188 


нччэоча OL NOILINOOULNI 


10856166500”- (9TINTIX LS 
089t£01820'- (SU NTIX 9S 
1141 95690"- (PI)NTIX SG 

У8199061”- (ET TÍNTIX PS 
8:260212"- (ZTINTIX €6 
8669091 °- (11) МТЇХ С6 
t0stoSbb^- (01) МЛ1Х 16 
1Є18620970- з (6)м1ЇХ 06 
1218626170- з (8)МТ1Х 6» 
608862071- = (L)NTIX 8p 
Б(99606'1- = (9) МТІХ L? 
6251559 °1- = (5) МТІХ 9р 
2962101 °2- = (р) МТІХ СЬ 
ErIEOOL'Z- = (EINTIX dp 
6118686"%- = (Z)NTIX Cb 
T9ETOPZ’S- = (TINTIX Zr 
990019667 = (91)1Х ТР 
16182216" 0% 
94518266" 6t 
22014418" 8t 
21806808" tE 
66800621" 96 
81108099" st 
92905196" = РЕ 
14411413 єє 
208616667 ze 
19166012" ТЕ 
88190161" 06 
8462221” 62 
Р%811490” 92 
692141400” tl 
$€S6672S00° = 92 
(ДУМАТУМ «| (91) уЗТҮНМ SZ 
(7) М31ЧН (СТ) Ұ3ТҰМ bz 
(С) М47ҰМ (ьт)матчн єс 
(y) VITVN (€) v4TVH ZZ 
(S) 4 ТҮН (ZU}WATWH TZ 
(9) ҮЗ ТҮН (11) У4ТҮН 02 
(4) Ч3ТҮН (01) ҮЗТҮН 61 
(8) уЗТҮН (6) ч 4ТҮН 81 

60862160" (8)ч АТЧН (1 
01:106160" (L)N ATWN 9I 
09281650980” (9) ч ITWN ST 
S66L6L* LO" (сум JTWN pt 
589Уу16290" (FIN JINN 
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Program "nbl6" (cont.) 
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Program "nbl6" (cont.) 
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